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The  removal  of  heat  from  fruits  and  vegetables  as 
soon  after  harvest  as  practical  has  long  been  recognized 
as  effective  in  retarding  the  ripening  process  and  con- 
trolling microbial  processes.   Presently,  very  little 
citrus  is  being  precooled.   Hydrocooling  of  citrus  has 
been  attempted,  but  was  not  satisfactory  due  to  the 
detrimental  effect  of  the  chemically  treated  water  on  the 
surface  appearance  of  the  fruit.   It  is  believed  that 
forced-air  precooling  of  citrus  would  receive  favorable 
attention  if  additional  basic  information  were  available 
on  heat  transfer  characteristics  of  citrus  in  bulk. 


XI 


Equations  describing  the  heat  transfer  in  a  bed 
of  citrus  fruit  were  derived.   A  numerical  solution  of 
these  equations  was  obtained  and  evaluated  utilizing  a 
digital  computer.   Temperature  was  determined  as  a 
function  of  1)  position  within  the  individual  fruit, 
2)  distance  from  inlet  of  bed,  and  3)  time.   Variable 
input  data  for  the  computer  program  include :   diameter 
of  fruit,  thermal  diffusivity  of  fruit,  specific  heat  of 
fruit,  density  of  fruit,  initial  fruit  temperature,  air 
temperature,  air  flow  rate,  total  weight  of  fruit,  depth 
of  bed,  and  convective  heat  transfer  coefficient.   An 
equation  was  developed  for  determining  the  average  con- 
vective heat  transfer  coefficient  for  a  bed  of  citrus  in 
the  form  of  a  Nusselt  number-Reynolds  number  relationship. 
The  constants  in  this  equation  were  experimentally  eval- 
uated for  oranges  and  grapefruit  utilizing  a  method  which 
made  use  of  the  computer  solution. 

Fifteen  experimental  cooling  tests  were  conducted 
on  oranges  and  grapefruit  in  bulk  containers.   Mass- 
average  and  surface  temperatures  were  measured  at  four 
levels  within  the  bed  of  citrus  with  air  velocity  ranging 
from  15  to  400  fpm. 

A  comparison  between  the  experimental  and  theoret- 
ical results  indicated  good  agreement  except  for  surface 
temperatures  of  grapefruit.   The  experimental  surface 
temperatures  of  grapefruit  were  consistently  lower  than 
those  predicted  by  the  numerical  solution.   This  was 
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shown  to  be  caused  primarily  by  the  low  thermal  conductiv- 
ity of  the  thick  rind.   A  solution  was  obtained  for  the 
temperature  gradient  in  a  composite  sphere  which  indicated 
a  much  steeper  temperature  gradient  through  the  rind  than 
through  the  juice  vesicles. 
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INTRODUCTION 

The  removal  of  field  heat  from  fruits  and  vegetables 
prior  to  processing,  transportation  and  storage  has  long 
been  recognized  as  effective  in  retarding  the  ripening 
process  and  controlling  microbial  processes  (21,  40,  41)—/. 
The  term  "precooling"  is  commonly  used  to  describe  this 
rapid  cooling  process  (1,  2,  42).   Maxie  et  al.   (34)  have 
described  the  biological  basis  for  prompt  cooling,  listing 
a  series  of  physiological,  biochemical,  pathological  and 
physical  changes  that  are  temperature-dependent  in  fruits. 
Precooling  is  likely  the  most  important  of  all  the  opera- 
tions used  in  establishing  and  maintaining  desirable,  fresh, 
and  salable  produce. 

According  to  Guillou  (22)  ,  the  rate  of  various  changes 
such  as  moisture  loss  and  the  action  of  enzymes  and  micro- 
organisms, are  commonly  doubled  or  tripled  with  each  10  F 
rise  in  temperature.   He  estimated  that  some  fresh  fruits 
and  vegetables  deteriorate  as  much  in  an  hour  at  90  F  as 
in  a  day  at  50  F  or  in  a  week  at  32  F. 

According  to  the  1972  Citrus  Summary  (17)  ,  approximately 
20%  of  Florida's  citrus  crop  is  consumed  in  the  fresh  form. 
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Numbers  in  parentheses  refer  to  numbered  references 


The  total  value  of  this  fruit  is  approximately  $110  million 
annually.   This  amounts  to  over  7%  of  the  total  value  of 
Florida's  annual  agricultural  production. 

At  present  most  of  the  citrus  fruits  are  not  cooled 
until  they  are  placed  in  the  refrigerated  vehicle  that  will 
carry  them  to  market.   The  refrigeration  systems  on  these 
vehicles  are  not  designed  to  remove  great  quantities  of  heat 
rapidly;  therefore,  the  temperature  of  the  fruit  is  not 
decreased  to  the  desired  level  for  a  considerable  period 
after  harvesting. 

It  has  been  conservatively  estimated  that  10%  of  all 
fresh  citrus  fruits  shipped  are  never  consumed  because  of 
waste  and  spoilage  (48)  .   The  value  of  these  fruits  is  over 
$10  million. 

A  small  quantity  of  citrus  fruit  is  presently  being  pre- 
cooled  by  hydrocooling  and  air  cooling.   Hydrocooling  of 
citrus  has  been  tried  by  some  producers,  primarily  due  to 
the  successful  hydrocooling  of  peaches  and  vegetables. 
However,  it  soon  became  apparent  that  fungicide  in  the  cool- 
ing water  was  essential  to  hold  decay  to  tolerable  levels. 
Unfortunately,  the  chemically  treated  water  had  a  detrimental 
effect  on  the  surface  appearance  of  the  fruit.   Also,  fruit 
cooled  with  water  seemed  more  susceptible  to  decay  upon 
warming  than  fruit  which  had  received  no  cooling.   Certain 
types  of  citrus,  notably  grapefruit,  exhibited  subsequent 
chilling  injury  (28) . 


It  is  believed  that  forced-air  precooling  of  citrus 
would  receive  more  favorable  attention  if  additional  basic 
information  were  available  on  heat  transfer  characteristics 
in  beds  of  citrus  fruit.   Much  of  the  research  on  precooling 
has  been  done  in  production  operations  or  in  experimental 
facilities  which  were  not  designed  for  obtaining  basic  heat 
transfer  data.   Thus  the  objective  of  this  investigation  is 
to  obtain  and  test  a  mathematical  model  that  can  predict 
the  temperature  distribution  within  a  bed  of  citrus  during 
forced-air  cooling.   Other  objectives  incident  to  this 
investigation  are  to  devise  methods  for  obtaining  thermal 
and  physical  properties  for  beds  of  citrus  fruit. 


REVIEW  OF  LITERATURE 
Studies  on  the  Thermal  Propertiesof  Fruits  and  Vegetables 

The  response  in  temperature  of  a  body  to  a  given  input 
or  removal  of  energy  in  the  form  of  heat  is  dependent  upon 
certain  properties  of  the  body.   These  properties  govern 
the  storage  and  transfer  of  heat  energy  through  the  body. 
Two  of  these  properties  are  specific  heat  capacity  (c) , 
BTU/lb-°F,  and  thermal  conductivity  (k) ,  BTU/hr-°F-f t . 
Another  property,  which  is  convenient  to  use  since  it  con- 
sists of  a  combination  of  three  properties  as  they  appear 
in  the  heat  conduction  equation,  is  thermal  diffusivity 
(k/pc),  ft2/hr. 

Since  most  fruits  and  vegetables  contain  a  high  per- 
centage of  water  it  can  be  expected  that  their  properties 
will  have  values  near  that  of  water.   Formulas  for  the 
specific  heat  and  thermal  conductivities  of  fruits  and 
vegetables  based  on  their  percent  water  have  been  developed 
(15). 

c  =  p/100  +  0.2  (100-p)/100  BTU/lb-°F  [l] 

k  =  0.32  /100  +  0.15  (100-p)/100  BTU/hr-ft-°F  [2] 


where 

c  =  specific  heat 

p  =  percent  water 

k  =  thermal  conductivity 

Gane  (20)  investigated  the  thermal  diffusivity  of  fruits 
and  vegetables  by  using  the  graphical  method  of  Gurney  and 
Lurie  (23)  for  solutions  to  the  conduction  heat  transfer  equa- 
tion for  various  shaped  bodies.   The  temperature  history  of 
of  the  center  of  fruits  and  vegetables  was  monitored  with 
thermocouple  junctions  during  cooling  and  was  used  in  con- 
junction with  these  graphical  solutions  to  estimate  the  ther- 
mal diffusivity.   He  found  the  following  average  diffusivity 
values  (ft2/hr) :   apples  --  0.0049,  oranges  —  0.0049, 
grapefruit  —  0.0049,  and  squash  —  0.0044.   He  also  found 
the  thermal  conductivity  of  oranges  and  grapefruit  to  be  0.24 
and  0.25  BTU/hr-f t-°F  respectively. 

Kethley  et  al.  (30)  conducted  a  series  of  cooling  experi- 
ments with  eight  different  fruits  and  vegetables.   They  used 
the  same  graphical  methods  of  Gurnie  and  Lurie.   In  the 
temperature  range  80  to  32  F  their  values  of  diffusivity 
ranged  from  0.00535  to  0.00615  ft2/hr. 

Turrell  and  Perry  (50)  found  the  specific  heat  by  using 
a  calorimeter.   They  obtained  values  of  0.885  for  Marsh 
grapefruit,  0.892  for  Eureka  lemons,  0.840  for  Valencia 
oranges  and  0.875  for  Washington  Navel  oranges. 

Several  investigators  have  determined  the  thermal  conduc- 
tivity (k)  of  citrus  by  assuming  it  to  be  homogeneous.   Turrell 


and  Perry  (50)  found  the  thermal  conductivity  averaged  0.23 
BTU/hr-f t-°F  for  Marsh  grapefruit,  Eureka  lemons,  Valencia 
oranges  and  Washington  Navel  oranges  when  all  varieties  were 
considered.   Perry  et  al.  (37)  determined  the  following  mean 
values  for  thermal  diffusivity:   Valencia  oranges  —  0.00368, 
Washington  Navel  oranges  —  0.00415,  Marsh  grapefruit  — 
0.00355,  and  Eureka  lemons  —  0.00417  ft2/hr. 

Smith  et  al.  (47)  developed  a  technique  of  evaluating 
thermal  diffusivity  which  corrects  for  deviation  from  the 
conventional  shapes.   The  technique  incorporates  a  geometry 
index  (G)  into  the  basic  Fourier  heat  conduction  equation  for 
a  sphere.   The  geometry  index  was  obtained  from  a  measure  of 
two  orthogonal  areas  of  the  shape.   Using  this  technique 
they  determined  an  average  thermal  diffusivity  of  0.00363 
ft  /hr  from  15  tests  run  on  Marsh  grapefruit.   Bennett  et  al. 

(7)  used  the  same  procedure  and  obtained  an  average  thermal 

2 
diffusivity  of  0.00355  ft  /hr  for  five  maturity  groups  of 

Marsh  grapefruit.   This  value  was  obtained  from  temperature 
response  at  the  approximate  mass  average  point.   Values  cal- 
culated from  temperature  response  at  the  center  and  one-half 
radius  point  were  slightly  lower,  probably  due  to  the  non- 
homogeneity  of  the  product. 

Bennett  et  al.  (7)  also  made  separate  determinations  for 
the  thermal  properties  of  the  rind  and  juice  vesicle  compo- 
nents of  Marsh  grapefruit.   Table  1  indicates  average  values 
for  each  component.   The  thermal  diffusivity  determined  here 


Table  1.   Comparison  of  Rind  and  Juice  Vesicle  Components 
of  Marsh  Grapefruit. 


Specific  Thermal       Thermal 

Heat        Density    Conductivity   Diffusivity 
BTU/lb-°F    lb/ft3     BTU/hr-ft-°F   xlO-3  ft2/hr 


Juice 

Vesicle 

0.907 

63.2 

0.2716 

4.56 

Rind 

0.829 

36.1 

0.1397 

4.67 

is  considerably  higher  than  that  determined  by  assuming  the 
fruit  to  be  homogeneous.   It  is  interesting  to  note  that  the 
thermal  diffusivity  of  the  rind  and  juice  vesicle  components 
has  values  very  near  the  same  while  the  thermal  conductivity 
and  density  variation  is  much  greater.   This  has  been  explained 
by  Perry  et  al.   (37)  that  in  a  given  material  where  density 
changes  because  of  changes  in  porosity,  the  thermal  conduc- 
tivity is  approximately  proportional  to  density,  so  that  the 
diffusivity  remains  about  constant. 

Bennett  et  al.   (7)  conducted  similar  tests  on  Valencia 
oranges  and  found  the  combined  thermal  conductivity  to  be 
0.248  BTU/hr-ft-°F. 

Studies  on  Fruit  and  Vegetable  Precooling 


Presently,  most  citrus  fruit  is  not  cooled  until  it  is 
placed  in  the  refrigerated  vehicle  that  will  carry  it  to 
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market.   The  refrigeration  systems  of  these  vehicles  are  not 
designed  to  remove  great  quantities  of  heat  rapidly;  thus 
the  temperature  of  the  fruit  is  not  decreased  to  the  desired 
level  for  several  days  after  harvesting.   During  this  time 
respiration  reduces  the  quality,  nutrient  value  and  shelf 
life  of  the  fruit,  thus  indicating  the  importance  of  pre- 
cooling. 

During  the  respiration  process  oxygen  from  the  air  is 
combined  with  carbon  from  the  plant  tissues  to  form  various 
decomposition  products  and  eventually  carbon  dioxide  and 
water.   In  addition  to  this,  the  enzymes  present  act  on 
various  substances  in  the  plant  tissues  and  gradually  cause 
changes  in  color,  texture,  and  chemical  composition  which 
mature  fruit  and  may  eventually  cause  serious  deterioration. 

During  the  respiration  process  heat  is  released  vary- 
ing in  amount  with  the  commodity  and  its  temperature. 
Table  2  illustrates  the  large  difference  between  the  heats 
of  respiration  at  high  and  low  temperatures. 

The  respiration  process  is  generally  believed  to  cause 
consumption  of  a  simple  hexose  sugar  glucose.   With  the  aid 
of  an  enzyme  system  and  oxygen,  glucose  sugar  is  reduced  to 
carbon  dioxide  and  water  in  accordance  with  the  following 
reaction: 


CrHn_Oc    +    60„    =    6C0_    +    6H_0   +    Heat  [3] 

6    12    o  2.  Z  2, 


Table  2.   Approximate  Heat  of  Respiration  of  Oranges  and 
Grapefruit  at  Various  Temperatures  (26) 


(BTU  per  24  hours  per  ton  of  fruit) 


Temp.  ( 

°F) 

Sensible 

Latent 

Total 

Heat 

Heat 

Heat 

Oranges 

32 

900 

63 

963 

40 

1400 

97 

1497 

60 

5000 

348 

5348 

80 

8000 

557 

8557 

Grapefruit 

32 

460 

32 

492 

40 

1070 

74 

1144 

60 

2770 

193 

2963 

80 

4180 

291 

4471 

The  enzyme  acts  as  a  catalyst  but  is  not  consumed  during 
the  reaction.   The  effectiveness  of  the  enzyme  as  a  catalyst 
is  greatly  reduced  at  low  temperatures  which  causes  a  reduc- 
tion in  the  heat  of  respiration  when  the  temperature  is 
reduced. 

Precooling  usually  refers  to  the  rapid  removal  of  heat 
from  a  fruit  or  vegetable.   In  regards  to  what  "pre"  means, 
Thevenot  (49)  stated:   a)  in  a  restrictive  sense:   previous 
to  transporation;  b)  in  a  wide  sense:  previous  to  transporta- 
tion or  to  storing;  c)  in  a  very  wide  sense:   previous  to 
transportation,  storing,  or  processing.   Redit  and  Hamer  (42) 
referred  to  precooling  as  the  rapid  cooling  of  a  commodity 
to  a  suitable  transit  or  storage  temperature  soon  after  har- 
vest, before  it  is  stored  or  moved  in  transit. 
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A  small  quantity  of  citrus  fruit  is  now  precooled  by 
several  methods.   Of  these  hydrocooling  and  air  cooling  are 
probably  the  only  two  that  can  be  combined  with  the  mechanized 
methods  of  packing  that  are  now  used.   However,  some  investi- 
gations on  hydrocooling  of  oranges  have  indicated  that  an 
increase  in  respiration  rate,  over  that  of  uncooled  oranges, 
occurred  when  they  were  allowed  to  return  to  room  temperature. 
Thus,  hydrocooling  can  actually  be  detrimental  to  oranges  (14) . 

On  the  other  hand,  Hopkins  and  Loucks  (28)  found  that 
air  precooled  oranges  were  outstanding  in  their  keeping 
quality  and  resistance  to  decay. 

Sainsbury  (43)  investigated  high  velocity  air  as  a  means 
of  rapidly  precooling  fruits.   In  cooling  experiments  involv- 
ing such  fruits  as  cherries,  apples,  and  apricots,  predicta- 
bility of  the  cooling  behavior  was  established  by  identifying 
the  cooling  rate  (CP)  for  a  given  experiment. 


Fruit  temperature  reduction 

ID.     -    TD 
i 

loge (TDi/TDf ) 


TD.  -  TDf 
CR  = Time       U] 


where 


TD .  =  initial  temperature  difference  between  fruit  and 

air. 
TDf  =  final  temperature  difference  between  fruit  and  air. 

This  equation  for  cooling  rate  (CR)  reduces  to 

TD. 
l09e-  TDT 

CR  "    Time  I5] 
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when  the  air  temperature  is  maintained  constant.   The  slope 
of  the  line  established  by  plotting  loge  (TDi/TDf)  as  a 
function  of  time  identifies  the  cooling  rate  (CR) . 

Guillou  (22)  reports  that  Newton's  law  of  cooling  will 
predict  quite  well  the  average  temperature  progression  of 
fruits  and  vegetables  during  cooling.   Newton's  law  is  appli- 
cable for  an  object  placed  in  surroundings  at  a  constant 
lower  temperature  provided  the  resistance  to  heat  transfer 
is  constant.   The  equation  expressing  this  law  using  the 
notation  of  Guillou  is 

dT  =  c  {T_T^  [6] 

where 

t  =  time  of  exposure  to  surroundings 

T  =  temperature  of  the  object  at  time  "t" 

T   =  temperature  of  the  surroundings 

C  =  "cooling  coefficient" 
The  integration  of  the  equation  by  separation  of  variables 
shows  that 

C  "  I  **.  Ipfe 

where  Tn  equals  the  initial  temperature  of  the  object.   This 
equation  for  the  "cooling  coefficient"  corresponds  to  the 
equation  of  Sainsbury  (43)  for  cooling  rate  (CR) . 

Bennett  et_al.  (  9  )  conducted  tests  on  forced-air  pre- 
cooling  of  Florida  citrus.   These  tests  were  conducted  on 
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citrus  in  bulk  containers  through  which  they  forced  high 
velocity  air.   They  developed  an  equation  to  predict  the 
cooling  rate.   Since  the  inlet  air  temperature  was  not  held 
constant,  the  air  temperature  at  the  end  of  the  run  was  used 
as  the  reference  temperature.   They  did  not  investigate  the 
temperature  distribution  within  the  bed. 

Leggett  and  Sutton  (33)  also  conducted  precooling  tests 
on  citrus  fruits.   They  developed  cooling  curves  and  obtained 
cooling  coefficients  for  oranges  and  grapefruit  in  crates. 
Both  hydrocooling  and  air  cooling  tests  were  made  for  dif- 
ferent crate  spacings  and  arrangements. 

Several  investigators  (5 ,  27 ,  36 ,  53 /   29)  have  con- 
ducted tests  involving  transport  phenomena  of  heat  and  mois- 
ture in  beds  of  biological  materials.   Most  of  these  investi- 
gations involved  analysis  similar  to  -chat  of  Schumann  (45) 
and  Furnas  (19 )  who  were  evidently  the  first  to  solve  the 
differential  equations  governing  heat  transfer  in  a  deep  bed 
of  particles.   For  a  system  in  which  a  fluid  passes  through 
a  prism  of  crushed  material,  Schumann  developed  a  series  of 
temperature  history  curves  subject  to  the  following  assump- 
tions : 

a)  that  the  particles  are  small  and  no  temperature 
gradient  exists  within  any  particle  at  any  time. 

b)  that  no  heat  is  produced  by  particles. 

c)  that  no  heat  is  transferred  from  particle  to  parti- 
cle by  conduction. 
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d)  that  the  fluid  is  incompressible. 

e)  and  that  the  thermal  properties  of  the  particles 
are  uniform  throughout  and  do  not  vary  with  time, 

Furnas  found  that  Schumann's  theoretical  curves  would  pre- 
dict the  temperature  behavior  of  beds  of  coke  and  iron 
through  which  air  was  moving. 


CONDUCTION  HEAT  TRANSFER  IN  SINGLE  CITRUS  FRUITS 

Theory  and  Assumptions 

Foundation  for  the  theory  of  heat  conduction  is  univer- 
sally attributed  to  Fourier  (18).   It  is  apparent  that 
Fourier  entertained  a  hope  that  practical  use  might  be  made 
of  his  work.   He  stated,  "It  is  easy  to  judge  how  much  these 
researches  concern  the  physical  sciences  and  civil  economy, 
and  what  may  be  their  influence  on  the  progress  of  the  arts 
which  require  the  employment  and  distribution  of  heat." 
Fourier's  work  was  awarded  a  generous  prize  by  the  French 
Royal  Academy  of  Science  in  1812.   Fourier's  work  was  treated 
like  a  work  of  art  for  more  than  100  years  with  no  practical 
use  being  made  of  it,  mainly  due  to  the  difficulty  encoun- 
tered in  obtaining  practical  solutions.   However,  in  1923 
engineers  Furney  and  Lurie  (23)  presented  their  solutions 
based  on  theory  but  with  an  experimental  character.   Answers 
to  many  practical  problems  of  interest  to  engineers  were 
made  readily  accessible  by  means  of  the  time-temperature 
charts  of  Gurney  and  Lurie.   Other  such  charts  have  been 
developed  by  means  of  the  experimental  technique  of  the  elec- 
tric analog  (25).   A  large  collection  of  time-temperature 
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charts  for  conduction  in  various  shapes  with  several  dif- 
ferent boundary  conditions  is  given  in  a  book  by  Schneider 
(44). 

A  considerable  amount  of  effort  has  been  given  by 
researchers  (24,  27,  38,  39,  46,  52)  to  correlate  the 
theoretical  solutions  for  regular-shaped  objects  with  exper- 
imental data  for  biological  materials  that  can  be  approxi- 
mated by  these  regular-shapes. 

Homogeneous  Sphere 

Most  citrus  fruits  have  the  general  shape  of  a  sphere 
and  therefore  could  be  expected  to  exhibit  heat  transfer 
characteristics  similar  to  that  of  a  sphere.   As  noted  in 
the  preceding   section,  citrus  fruits  are  not  composed  of 
homogeneous  materials  but  consist  of  two  main  sections,  the 
rind  and  the  juice  vesicles.   It  was  also  noted  that  although 
the  thermal  conductivity  (k)  and  density  (p)  were  quite 
different  for  the  two  portions,  their  thermal  dif fusivities 
(— )  were  surprisingly  close  to  one  another. 

Bennett  et  al.  (8)  investigated  the  temperature  distri- 
bution in  Marsh  grapefruit,  which  has  a  rind  that  constitutes 
30  to  40  percent  of  its  total  volume,  and  found  that  the 
presence  of  the  rind  has  very  little  effect  on  the  tempera- 
ture profile.   This  is  indicated  in  Figure  1  which  shows  the 
internal  temperature  distribution  in  Marsh  grapefruit  plotted 
from  raw  experimental  data  recorded  during  transient  cooling. 
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on 


•0.4        C.j 

RADIUS  RATIO 


Figure  1.   Internal  temperature  distribution  in  Marsh 

grapefruit  plotted  from  raw  experimental  data 
recorded  during  transient  cooling  test  runs. 
Dotted  lines  show  discontinuity  in  temperature 
gradient  at  interface. 
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The  dotted  lines  show  a  discontinuity  in  the  temperature 
gradient  at  the  interface  between  the  juice  vesicles  and 
rind. 

On  the  basis  of  these  findings,  it  will  be  assumed  that 
citrus  fruits  may  be  treated  as  homogeneous  without  signifi- 
cant error  in  the  temperature  profile. 

Heat  of  Respiration 

The  justification  for  not  including  the  heat  of  respira- 
tion as  a  heat  source  in  the  conduction  equation  is  bj.sed  on 
the  findings  of  Awberry  (4)  and  Hood  (27).   Awberry  calculated 
the  temperature  excess  at  the  center  of  an  apple  due  to  its 
heat  of  respiration,  after  two  hours  of  cooling  and  found 
it  to  be  0.0232  C,  while  Hood  calculated  a  value  of  0.011  F 
at  the  center  of  the  cucumber  after  30  minutes  of  cooling. 
Both  of  these  products  have  heats  of  respiration  higher  than 
citrus,  thus  the  effect  of  heat  of  respiration  on  citrus 
during  cooling  is  believed  to  be  negligible.   It  should  be 
pointed  out  that  the  effect  of  heat  of  respiration  should 
not  be  neglected  in  some  other  cases  such  as  storage  at 
elevated  temperatures. 

Analytical  Solution 

The  governing  equations  for  conduction  heat  transfer 
from  a  homogeneous  sphere  without  heat  sources,  initially 


at  a  uniform  temperature  (T,  )  and  surrounded  by  a  fluid  of 
constant  temperature  (T  )  are 


3T     ,32T    2  3T,  r  Q, 

3t  =  a(872  +  r   3^  [  8] 


T(r,o)  =  T,  [9] 


(k  i)r=a  =  h(Ts  "  T~>  [10] 


where 

t  =  time,  hr 

r  =  radial  distance  from  center,  ft 

a  =  radius  of  sphere,  ft 

T  =  temperature  at  r  and  t,  °F 

2 

a  =  thermal  diffusivity,  ft  /hr 

k  -   thermal  conductivity,  BTU/hr-ft-°F 

The  theory  for  the  following  solution  of  equations 

[8],  [9],  and  [10]  was  taken  from  Arpaci  (  3)  and  equation 

[21]  agrees  with  his  solution. 

Equation  [  8 ]  can  be  written  as : 


3T     a   3  ,  2  3T 


and  substituting  e  =  T  -  T   into  equations  [9] ,  [10]  ,  and 

oo 

[11]  gives: 


39  -       3,236,  r.0, 

at  "  ?2-  37  (r  5r>  [12] 
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(r,0)  =  q1  [13] 


-k  [56,(^t}  ]  =  hG(a,t)  [14] 

a  r 


One  more  boundary  condition  is: 

6(0, t)  =  finite  or  dd^'t]    =  0        [15] 
Using  the  transformation: 

6  (r,t)  =  y  (r,t)/r  [16] 

reduces  the  spherical  Laplacian  to  the  cartesian  Laplacian 
whose  solution  is  expressible  in  terms  of  circular  functions 

We  may  now  express  equations  [12]  through  [15]  in 
terms  of  ¥  by  using  equation  [16].   The  result  is: 


ff-a^,       fCr,0)-r 

3t      3r 


[17] 


Y(0,t)  =  0,       -k  3^'t}  =  (h  -  |)  <F(a,t) 


Hence  the  problem  is  reduced  to  a  problem  of  cartesian  geome- 
try. 

The  product  solution  ^(r/t)  =  R(r)  t  (t)  yields: 


^|+  A2  R=  0 
dr 


R(0)  -  0,     d?(a)  +  (^  -  -)  R  (a)  =  0      [18] 
'     ar      k   a 
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in  r  and  equation  [19]  in  time. 

JL  +  aA2   =0  [19] 

dt 

The  solution  of  equation  [18]  is: 


R  (r)  =  A*  (r),     4>  (r)  =  sinA  r 
n       n  n         n  n 


and  the  zeros  of 


(A  a)  cosA  a  =  (1-Bi)  sinA  a 
n       n  n 


are  the  characteristic  values,  where  Bi  =  ha/k, 
The  solution  of  equation  [19]  is: 


x   t)  =  C  e   n 
n       n 


Thus  the  product  solution  becomes: 

A  2<- 
¥(r,t)  =  I  be  a  n  fc  sinA  r       [20] 

n  n 

n=l 
The  initial  value  of  equation  [20]  is: 


r0,  =  Z     b   sinA  r 
1         n      n 

n=l 

The  coefficient  b   is  obtained  by  the  expansion  of  r6 ,  into 

n 

a  Fourier  sine  series  evaluated  from  zero  to  a. 
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20,  (sinA  a  -  X  a  cosA  a 

1  n n n 

n    A   (A  a  -  sinA  a  cosA  a) 
n    n        n      n 


Thus  the  unsteady  temperature  of  the  sphere  is: 


00  2 

T(r,t)  -  T  sinA  a  -  A  a  cosA  a    -aA  t  sinA  r 

T,  -  T   "  *X  a  -  sinA  a  cosA  a;  e        A  r 

1    °°         ,    n        n      n  n 

n=i 

[21] 


Numerical  Solution  for  Heat  Conduction  in  a 
Sphere  with  Time-Dependent  Boundary  Conditions 


The  analytical  solutions  for  the  heat  conduction  equa- 
tion presented  in  the  preceding  section  are  based  upon 
assumed  initial  and  boundary  conditions  which  can  seldom  be 
satisfied  in  practical  applications.   Especially  in  dealing 
with  many  biological  products,  there  is  a  need  for  informa- 
tion when  initial  and  boundary  conditions  may  exist  such  as: 
with  time-varying  environmental  conditions  or  surface  con- 
ductances, with  change-of-phase  or  mass  transfer  at  the 
boundary,  or  with  initial  object  temperature  distributions. 

The  only  practical  method  for  solving  this  type  problem 
is  by  numerical  methods,  such  as  those  presented  by  Arpaci  (3) 
Since  cooling  in  a  deep  bed  involves  time-dependent  boundary 
conditions,  a  numerical  solution  of  the  heat  conduction 
equation  for  a  sphere  surrounded  by  a  fluid  with  arbitrarily 
varying  temperature  will  be  derived. 
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The  difference  equations  which  approximate  the  govern- 
ing differential  equations  can  be  obtained  by  approximating 
each  term  of  the  differential  equation  with  its  equivalent 
finite  difference  expression  or  by  performing  a  heat  balance 
method  on  each  type  of  node  in  the  solid.   The  heat  balance 
method  will  be  used  in  this  formulation.   For  a  sphere  there 
are  three  types  of  nodes  as  indicated  in  Figure  2,  surface 
nodes,  interior  nodes,  and  center  nodes. 

The  first  law  of  thermodynamics  applied  to  the  interior 
node  gives: 

qi_1  [4?r  (r-^f)  2]+qi+1  [4ti  (r+^f )  2]  = 

Tn+l_Tn 

_  4  ,.  ,  Ar.  3  ,   Ar,3,   i  i      r_.0l 

pC  ^[(rH—j)  -  (r — ~)     ]  — ^ [22] 

where 

T    denotes  the  temperature  at  the  ith  node  at  time  n+1, 
Now  stating  Fourier's  law  of  conduction,  we  have: 

Tn    ,-Tn  T?+1-T? 

.,     i-l      i  i       i+l      i 

qi+l   =   K  Ar"^    '         qi+l   =   k  Ar~ 

Substituting  into  equation    [22]    results    in: 

Tn   _„n  Tn   _   n 

4lT^  f-r-  HL\2      i-i   i    ilk  ,   Ar,  2   i+l    i  . 
"Ar  <r~    2'  Ar       Ar  K,        2J  Ar 

Tn+l_Tn 

4   r.  ,Ar,3    ,   Ar.3,   i i 

pc  3  Tr[(r+-2-)   -  U~)     ]  ir- 
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Figure  2.   Diagram  indicating  a  heat  balance  on  the 
three  types  of  nodes  used  in  a  sphere 
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rearranging: 


Tn+1  m    {1  _  oAt   [(1_A£}2  +  (1+A|,2]}   Tn 


+  aA  t   rM_Ar,2  mn    ^    n  _,_ Ar  ,  2 


(Ar) 


2-  [d-^   Tx-1+  (1+iTr)   Ti+1]     [23] 


where   {    }  =  0   to  insure  convergence . 


Now  consider  the  surface  node.   The  first  law  of  thermo- 
dynamics applied  to  the  surface  node  gives: 

q._1[4,(R-A|)2]  +  qconv[4TTR2]  = 

Tn+1-Tn 
PC  !*[R3-(R~§>3J  -i-_L 

Substitution  for: 


Tn  .-Tn 

q.  .  =  k  -^ ,    q      =  h(Tn-Tn) 

^i-l        Ar       Hconv      °°   i 


gives : 


k  ^^i  [4U(R-Af)2]  +  h^-T^tilRih] 


4   rr,3  ,„  Ar,Jn 
pc  o-  7T  [R  -  (R T)  ] 


Tn+l_Tn 


3  " L"   x"   2'  J     At 


rearranging: 
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Ar  >  2 
Tn+1  =  3ctAt     U   2}       Tn    + 
i       Ar    [R3_(R_Ar)3]    i-1 


2 

3R  haAt n 

k[R3_(R_Ar)3]    . 


{1 3aA_t R"  2     +  S_ll  }    Tn    [24] 

[R3_(R_A|)3]      Ar 


where   {   }  >0  to  insure  convergence. 

Applying  the  first  law  to  the  center  node  results  in 

n+1  n 
,Arx2,       4   ,Ar.  3  ±i  i 


qi_1[47r  (—^)    ]    =    pc  —  -n  (— ~) 


At 


substituting  for: 

qi+l  =  k 


l+l   i 


gives 


.  Ti+l"Ti  r.  ,Ar,2,     „  4   ,Ar,3  T1?"1"1-^ 

k  — —  [4tt  (-^)  ]  =  pc  T  it  (-~)   _i i 

Ar        z  J     ^      At 


Simplifying  and  rearranging: 

Tn+1  =  6aAt    Tn    +   {l  _  .^oAt  }  T^      j25] 
1      (Arr  (Ar)'1 

where    {   }  >  0  to  insure  convergence. 

No  similar  solution  was  found  in  the  literature  for  comparison, 
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A  computer  program  was  developed  to  evaluate  equations 
[23]/  [24],  and  [25]  for  temperature,  T,  as  a  function  of 
"r"  and  "t."   The  arbitrarily  varying  fluid  temperature, 
T  ,  is  entered  as  input  data.   A  flow  diagram  for  the  com- 
puter program  is  given  in  Figure  3  and  the  computer  program 
is  shown  in  Appendix  A. 

In  order  to  test  the  accuracy  of  the  numerical  solution, 
the  air  temperature,  T  ,  was  entered  as  a  constant  so  that 
the  results  could  be  compared  with  the  exact  solution 
obtained  in  the  preceding   section.   This  comparison  is 
shown  in  Figure  4.   The  data  used  to  generate  these  curves 
are  the  following:   radius  =  0.113  ft,  k  =  0.256  BTU/hr-ft- 
°F,  h  =  11.7  BTU/hr-ft2-°F,  a  =  0.0051  ft2/hr,  c  =  0.9 
BTU/lb-°F. 
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Read:   a,  k,  h,  radius,  number  of  node  points  in  sphere, 
initial  temperature  of  fruit,  time  increment, 
total  cooling  time,  print  frequency,  test  number, 
arbitrarily  varying  air  temperature. 


1 

[Compute  combinations  of  input  data  which  remain  constant.! 

i 

I  Set  initial  conditions. 

' 

] Compute  air  temperature  from  input  data. 

1 

' 

Compute  surface  temperature  (Equation  [23]). 

i 

j Compute  interior  node  temperatures,  except  center 
(Equation  [24] ) . 

i 

Compute  center  temperature  (Equation  [25]) 

j  increase 
]  time  by  one 
increment 


Yes 


j Print  time, [ 
lair  temp, 
ifruit  temp 


Yes 


Figure  3.  Flow  diagram  for  digital  computer  program  of 
heat  conduction  in  a  sphere  with  arbitrarily 
varying  air  temperature. 
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O  Analytical  solution 
—  Numerical  solution 


-  —  Surrounding  air  temp 


CENTER 

MASS    AVERAGE 

SURFACE 


0-5  1.0 

COOLING  TIME-HOURS 


Figure  4.   Comparison  of  analytical  solution  with  numerical 
solution  for  heat  conduction  in  a  sphere 


TEMPERATURE  DISTRIBUTION  WITHIN  A  BED 
OF  CITRUS  DURING  COOLING 


Laplace  Transform  Solution  of  Simplified  Model 

In  the  preceding   section,  solutions  for  temperature 
distribution  within  a  single  fruit  were  obtained.   Now 
the  analysis  will  be  extended  to  include  temperature  dis- 
tribution within  a  bed  of  citrus  fruit.   In  this  section 
a  simplified  model  will  be  assumed  in  order  that  an  exact 
mathematical  solution  can  be  obtained.   This  simplified 
model  was  developed  by  making  the  following  assumptions: 

a)  No  mass  transfer  occurs  between  the  citrus  fruit 
and  the  cooling  medium. 

b)  There  is  no  temperature  gradient  within  the 
individual  fruit. 

c)  The  air  velocity  in  all  interstices  of  the  bed 
will  be  constant. 

d)  No  conduction  heat  transfer  between  the  fruit. 

e)  Thermal  properties  of  air  and  fruit  are  considered 
constant. 

These  assumptions  are  similar  to  those  made  by 
Schumann  (45) •   The  first  assumption  is  believed  to  be 
reasonable  due  to  the  small  amount  of  moisture  loss  from 
citrus  during  cooling  (less  than  1%  for  normal  conditions). 
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This  assumption  has  apparently  introduced  a  considerable 
error  for  studies  involving  fruits  and  vegetables  which 
do  not  have  a  waxy  surface. 

The  second  assumption  is  believed  to  be  the  weakest 
and  probably  causes  a  considerable  error  for  the  faster 
cooling  rates.   However,  this  error  is  believed  to  be 
minimized  by  using  the  mass-average  temperature  for  the 
product.   Although  the  last  three  assumptions  are  not  pre- 
cise, their  effect  is  believed  to  be  small. 

In  the  next  section  the  second  assumption  will  be 
investigated  by  obtaining  a  numerical  solution  which 
accounts  for  temperature  distribution  within  the  fruits. 
The  numerical  solution  will  be  compared  with  the  exact 
mathematical  solution  obtained  in  this  section,  for  a 
slow  cooling  rate  which  corresponds  to  a  negligible  temper- 
ature gradient  in  the  individual  fruit. 

In  order  to  obtain  the  differential  equations  describ- 
ing the  product  temperture  in  a  bed  of  citrus  fruit  and 
the  temperature  of  the  air  flowing  through  it,  a  heat 
balance  will  be  written  for  a  volume  A  dx  as  shown  in 
Figure  5.   "A"  is  the  cross-sectional  area  of  the  bed  and 
dx  is  the  height  of  the  control  volume.   The  following 
symbols  will  be  used  in  this  derivation: 

T  ,  temperature  of  air,  °F 
a 

T  ,  temperature  of  product,  °F 
h,  convective  heat  transfer  coefficient,  BTU/hr  ft  °F 
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Figure  5.   Control  volume  for  heat  balance  on  bed 
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A*,  product  surface  area  per  cubic  foot,  ft 

2 

A,  cross-sectional  area  of  the  packed  bed,  ft 

V  ,  superficial  velocity  of  air,  ft/hr 
a 

p  ,  product  bulk  density,  lb/ft 

3 

p  ,  air  density,  lb/ft 
a 

c  ,  specific  heat  air,  BTU/lb  °F 

c  ,  specific  product,  BTU/lb  °F 
P 

6  ,  dimensionless  temperature  ratio  of  the  air 

6  ,  dimensionless  temperature  ratio  of  the  product 

Y,  dimensionless  height  of  packed  bed 

Z,  dimensionless  time 

t,  time,  hr 

u,  absolute  viscosity,  lb/ft-hr 

x,  distance  from  air  inlet,  ft 
The  amount  of  air  flowing  past  a  point  "x"  in  the  amount  of 
time  dt  can  be  expressed  as: 


V  Apa  dt  (lb)  [26] 

a    a 


The  heat  content  of  the  air  entering  A  dx  in  dt  is 


V   Ap   c   T   dt  (BTU)  [27] 

a   a   a   a 


The  heat  content  of  the  air  leaving  A  dx  in  dt  is : 

3T 

V   Ap   c   (T   +  v-^  dx)  dt  (BTU)  [23] 

a    a   a    a    ox 
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The  change  in  enthalpy  of  the  air  with  respect  to  dx 
is  the  difference  between,  equations  [2  7]  and  [28]  or 


T 
Va  APa  ca  (  af  dx)  dt  (BTU)         ^29^ 


The  change  in  enthalpy  of  the  air  in  the  control  volume 
with  respect  to  time  is  equal  to: 


3T 
A  dx  Pa  ca  if  dt  t3°] 


The  total  change  in  enthalpy  of  the  air  is 


3T    3  T 
A  dx  pa  ca  dt  (Va-1|+  ^-f)        [31] 


As  air  moves  through  the  bed,  convective  heat  transfer 
occurs.   The  enthalpy  change  of  the  air  in  time  dt  that 
results  from  this  heat  transfer  is: 


qconv  =  h  A*  (Tp~Ta}  A  dx  dt        [32] 


where  A*  is  the  surface  area  of  product  per  unit  volume  of 
the  bed. 

Setting  equation  [31]  equal  to  [32],  the  enthalpy 
change  in  the  air  with  respect  to  time  and  position,  gives: 

8T       3T 
A  Pa  Ca  (lf  +  Va  ~&    dx  dt  = 


h  A*  (T   -T  )  A  dx  dt  [33] 

p    a 

or 
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'  +  --"  '  " (T   -  T  )  =  0   [34] 


V     St    3x.   V   p^  c      a     p 
a  a  "a   a 


Equation  [34]  is  the  result  of  writing  a  heat  balance  on  the 
air  in  the  control  volume  A  dx. 

Now  a  heat  balance  will  be  written  for  the  citrus  fruit, 
At  any  time  t  the  enthalpy  of  the  volume  A  dx  of  the  bed  is: 


p   A  dx  c   T  [35] 

P       P   P 


At  time  t  +  dt  the  enthalpy  is : 

Pp  Adx  cp  (Tp  +  ^  dt)  [36] 

The  change  in  enthalpy  of  the  product  in  time  dt  is 
obtained  by  substracting  equation  [36]  from  [35]: 

ST 

-p   A  dx  c  -£     dt  [37] 

p         p   at 

This  change  in  enthalpy  of  the  product  in  the  time  dt  is  the 
result  of  heat  transferred  from  the  product  (equation  [ 37] ) . 
Equating  equations  [37]  and  [33]  gives: 


-pp   CP  -st  -   h  A*  (TP  -  V        [38] 


Schumann  (45)  was  the  first  person  to  obtain  a  solution  of 
equations  [34]  and  [38].   Businger  (10),  Bakker  -  Arkema  and 
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and  Bickert  (5),  and  Huong  et  al .  (  29)  observed  that  the 

3T 

— ;r£-  in  equation  [34]  is  very  small  compared  to  V  ,  the  ve- 

3T  a 

locity  of  the  cooling  air,  thus  — ~     may  be  omitted  from 

a  t 

equation  [34]/  making  it  possible  to  obtain  a  solution  by 
Laplace  transform.   Then  the  two  equations  which  must  be 
solved  are: 


7)T 

a  h  A*     ,m    m  .  .   , 


a  "a   a 


and 


8Tp  +   h  A* 

3t    p   c    VJ"p    xa 
P   P 


(T   -  T  )  =  0  [40] 


The  theory  for  this  solution  and  the  Laplace  transforms 
were  taken  from  Churchill  (12). 

It  is  advantageous  to  make  the  variables  of  temperature, 
time,  and  position  dimensioniess  as  follows : 

Temperatures : 


T      -     (T    ) 

Tt~T -     (T    ) =    0a  C41] 


T      -     (T    ) 

p       up;t=o  _  u   , 

UaJx=0         {   p;t=0  p 
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Distance : 


h  A*  r 

x  =  Y  [43] 


p   c   V 
ra  a   a 


Time 


h  A* 

t  =  Z  [44] 


P   c 

P  P 


Substituting  equations  [41],  [42],  [43],  and  [44] 
into  equations  [39]  and  [40]  gives: 


30 

"37  +  Ga  "  0p  =  °  I  ^ 


30 

-r£  +0    -  0   =  0 

3Z     p     a 


[46] 


In  this  study  the  initial  fruit  temperatures  were  uniform 
and  the  cooling  air  temperature  was  constant.   This  is 
expressed  in  terms  of  boundary  and  initial  conditions  as: 


TP  <X'  0)  "  (Vt=o  [47] 


and 


Ta  (0,  t)  =  (Ta)x=Q  [48] 
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In  dimensionless  form  equations  [47]  and  [48]  becomes 


(T  )     -  (T  ). 

0   (Y,  0)  =    p  -°  _  ,?V-°  =  0     [49] 
y  v  ayx=o    v  p't=o 


and 


(Ta}x-0  -  (Vt-o 

©    (0,  Z)  =  -J^N^ /mP^  [50] 

a  (T  )     -  (T  ). 

a  x=o      p  t=o 


Equations  [45],  [46],  [49],  and  [50]  can  be  solved  by 
Laplace  transforms.   Let  Laplace  transforms 


Ly  [0a  (Y,  Z)]  =  F  (s,  Z) 


L   [0   (Y,  Z)]  =  G  (s,  Z) 

ir  1 

and 


Ly  [0a  (Y,  Z)]  =  SF  (s,  Z)  -  0a  (0,  Z) 


Taking  the  Laplace  transforms  of  equations  [45]  and  [46] 
w.  r.  t.  Y  (Y  -y   s)  gives: 

s  F  (s,  Z)  -  0   (0,  Z)  +  F  (s,  Z)  -  G  (s,  Z)  =  0 

cL 

[51] 

and 

dG  (s,Z) 


d  Z 


+  G  (s,  Z)  -  F  (s,  Z)  =  0     [52] 


38 


at 


Y  =  0,  03  (0,  Z)  =  1 
a 


substituting  in  equation  [51]  gives 


F  (s,  Z)  =  1  VVl'  2)  E53] 


substitute  equation  [53]  into  equation  [52],  then 


d  G  (s,  Z)      1      G  (s,  Z)       .     . 

'  d  z      i~TT  +  s  +  l G  (S/  Z)   ~ 


G  (s,  Z) 


or 


s  +  1    s  +  1 


d  G  (s,Z)      s       .     .   .    1        rul 

d~z +  i~TT  G  (s'  z)  ~  i^TT     [54] 


Equation  [54]  is  an  ordinary  differential  equation  of  first 
order;  the  general  solution  is: 


G  (s,  Z)  =  i  +  Cxe  -  s  ^  1   Z  [55] 


Using  I.C.  at  Z  =  0 


0,     (Y,  o)  =  0 

d 


Thus , 


39 
G  (s,  0)  -  0 
substitution  into  equation  [55] ,  yields 


1      s 


Gl,'Z)"?"ie"  s-f-1   Z       1563 


Taking  the  inverse  transform  w.  r.  t.  s.  (s  ->  Y)  ,  gives 


0p  (Y,  Z)  =  Ls_1  [G  (s,  Z)3  = 


Ls_1  t|  -  J  e  -  Hi z]  [57] 


where  Lg    is  the  symbol  for  the  inverse  Laplace  transform 

w.  r.  t.  s. 

11  s 

Now  rearrange    [-  -  —  e   -        +        Z]    so   that   available 

transform  pairs    can   be   used;    since 


1  e §_.     Z 

s  s   +   1 


1        1  1 „        1 

—  +   —  e    -    ; — r-   Z    -    — 

s  s  s    +    1  S 


ss   +   lls         }Le        s   +   lZ        1J 
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S   Z 


s    s  +  1   Q 


s  +  1 
/   e         dZ 


1       /•      1  S+1       -,  „ 


1       ,  1        S  +  1     -,-, 

a  "  /  e    bTT  e      dz 


Thus, 


s  1 

Z      Z    -Z  -tv  z 


1    1      s+1       ,  1      s+1 

—  -  —  e         =   /   e     ; — -  e       dZ 

s    s  0         s  +  1 


and  we  have : 


_      Z    -Z         -~  Z 
0(Y,Z)  =  L  "  [  /   e     — ~T  e  S+±        dZ  ] 

P  &      Q  i>  -r  x 


[58] 


Utilizing  the  substitution  theorem  and  tables  of  Laplace 
transform  pairs  given  by  Churchill  (12) : 


f(s-a)  -  L  [eat  F(t)]  [59] 


tv  c  k/s  -  L  [(t/k)  (V  1)/2  Iv_1  (2/kt)        [60] 


where  v>0 
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For  this  application,  a  =  1,  t  =  Y,  v  =  1,    k  =  Z.   Applying 
equation  [59]  to  equation  [58]  gives: 


0   (Y,  Z)  eY  =  L"1[(    /   e"Z  dZ) (    /  f   eZ/s  dZ)] 


Now  applying  equation  [60]  gives: 


0p  (Y,  Z)  =  e~Y      /   e"Z  Iq  (2  JTZ)    dZ   [61] 


which  is  the  solution  for  the  product  temperature  as  a 
function  of  position  in  the  bed  and  time. 

Differentiating  equation  [61]  w.r.t.  Z  gives: 

38P"'Z)   -a      /e-<Y+  Z'  1(2  ^X)  dz 

0  4  d  4  O 


=  e"(Y  +  Z)  IQ(2  /£x)  dZ 
substituting  into  equation  [46]  yields: 


0a(Y,  Z)  =  ~ g^'  Z)   +  0a  (Y,  Z) 

Z 
=  e-(Y  +  Z)  x     ,       y^)    +  e-Y     ;   e-Z  j  (2  ^2}  dz 
o  o 

0 

[62] 
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which  is  the  solution  for  the  air  temperature  as  a  function 
of  position  in  the  bed  and  time. 

In  equation  [61]  and  [62]  ,  I  (2  /Yz")  is  known  as  the 
modified  Be&sel  function  of  the  first  kind,  of  order  zero, 
with  argument  2  /YZ";  in  general, 


(x/2) 


Iv(x)  "  J-        k:  (k+v)  : 

K=o 


Equations  [61]  and  [62]  have  been  evaluated  numerically 
(19)  and  are  shown  in  Figures  6  and  7. 

Numerical  Solution  of  Realistic  Model 

The  notations  and  assumptions  used  in  the  preceding 
section  will  apply  here,  except  that  temperature  gradients 
within  the  individual  fruit  may  exist,  but  are  a  function  of 
radius  only.   Thus,  an  additional  equation  is  required  to  de- 
scribe the  conduction  heat  transfer  within  the  fruit.   This 
equation,  with  its  initial  and  boundary  conditions,  was 
discussed  in  a  previous  section  (equation  [3],  [9],  [10]). 

Consider  an  energy  balance  written  on  the  air  flowing 
through  a  differential  volume  a  dx,  as  depicted  in  Figure  5. 

Energy  out  =  energy  in  +  change  in  internal  energy  of 
the  product: 

3T 

VpcAtdt   +  VpcA  —A  dx   dt   = 
a      a      a  a  a    "a      a  ox 

3T 

VpcAtdt-pc       -7r—       ,.     ,.      , 
a    Ka      a  a  Kp      p      3r.      at  A   ax 


Q. 
0 


Figure  6.   Schumann's  curves  for  temperature  history 
of  solid  for  values  of  Y  from  0  to  10 


CD 


Figure  7.   Schumann's  curves  for  temperature  history 
of  gas  for  values  of  Y  from  1  to  10 
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or 


W.jE^L_!!g  [63] 

3x   V   p   c    at 
a   a   a 


Thus,  equations  [8],  i  9],  [10],  [6  3]  describe  the  heat 
transfer  within  a  bed  of  fruit  in  which  temperature  gra- 
dients exist  within  the  individual  fruit.   A  numerical 
solution  for  this  set  of  equations  will  be  developed 
utilizing  the  numerical  solution  that  was  previously 
obtained  for  the  conduction  equation. 

Considering  the  differential  layer  of  the  bed  to  be  a 
finite  layer,  equation  [63]  becomes: 


AT  p   c   AT 

a  MP   o    p                       r  c  /n 
=  — <■ « [6  4] 

Ax  Va  pa  ca  At 


Equation  [64]  along  with  equations  [23],  124],  [25]/  which 
were  derived  in  a  previous  section,  were  used  to  develop  a 
computer  program  which  results  in  a  printout  of  fruit 
temperature  as  a  function  of  1)  position  in  bed,  2)  point 
in  the  individual  fruit  and  3)  time.   In  order  to  solve 
this  set  of  equations  it  is  necessary  to  compute  the  mass- 
average  temperature  of  the  fruit  after  each  time  period. 
This  is  accomplished  by  assuming  a  linear  temperature  pro- 
file within  each  segment  of  the  fruit.   This  causes  the 
mass- average  temperature  to  be  located  at  the  mass-average 
point,  which  is  located  at  r/radius  equal  to  0.7937  for  a 
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sphere.   Mass- average  temperature  denotes  a  single  value 
from  the  temperature  profile  of  a  product  that  would  become 
the  uniform  product  temperature  upon  reaching  equilibrium 
under  adiabatic  conditions. 

Since  the  temperature  profile  in  a  sphere  is  not 
linear,  the  mass- average  temperature  will  not  be  located 
exactly  at  the  mass-average  point  and  will  shift  according 
to  the  shape  of  the  temperature  profile.   For  this  solution 
the  mass-average  temperature  of  the  fruit  was  determined  by 
computing  the  mass-average  temperature  of  each  spherical 
shell  and  then  computing  a  weighted  average  temperature  based 
on  the  mass  of  each  spherical  shell. 

In  order  to  account  for  the  shift  in  the  location  of 
the  mass- average  temperature  due  to  the  low  heat  capacity 
of  the  rind  (pC) ,  the  outer  and  inner  shells  of  the  spheri- 
cal model  were  given  values  for  heat  capacity  corresponding 
to  that  of  the  rind  and  juice  vesicles  as  determined 
separately.   The  heat  capacities  for  the  rind  and  juice 
vesicles  of  oranges  and  grapefruit  were  obtained  from  the 
findings  of  Bennett  et  al .  (7)  as  23.7,  57.6,  30.0,  and 
57.2  BTU/ft3-°F  respectively. 

The  average  thickness  of  the  rind  for  Valencia  oranges 
and  Marsh  grapefruit  was  also  taken  from  this  reference  and 
roughly  corresponds  to  the  outer  segment  of  a  ten-segment 
sphere  for  oranges  and  the  outer  two  segments  for  a  ten- 
segment  sphere  for  grapefruit. 
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Soule  et  al.  (4  8)  experimentally  located  the  points 
of  mass-average  teraperature  during  fcrced-air  cooling  for 
oranges,  grapefruit  and  tangeios.   The  values  of  r/radius 
were  0.7856,  0.7722,  and  0.7353,  respectively.   Evidently 
the  lower  densities  of  the  rind  shift  the  location  of  the 
mass  average  temperature  towards  the  center  of  the  fruit. 
Soule  et  al .  also  presented  their  results  in  terms  of  locus 
of  mass  average  temperature  versus  cooling  time  for  a 
specific  cooling  rate  (Figure  10)  .   Bennett  et  al .  (8) 
conducted  the  same  type  of  tests  for  grapefruit  cooled  in 
agitated  ice  water  (Figure  9  )  and  got  a  much  larger  shift 
in  the  location  of  the  mass-average  temperature. 

Each  step  of  the  numerical  solution  presented  in  this 
section  is  indicated  by  a  flow  diagram  for  the  computer 
program  (Figure  8).   The  computer  program  and  sample  output 
data  for  cooling  a  deep  bed  of  citrus  fruit  are  shown  in 
Appendix  B.   This  output  is  shown  in  graphical  form  in  the 
Data  and  Results  section  of  this  chapter. 

Determination  of  Ccnvective  Heat  Transfer 
Coefficient  and  Thermal  Diffusivity 

Theory 

Transfer  of  heat  by  transport  of  fluid  material  is  known 
as  convection.   The  transport  may  be  natural  convection, 
caused  by  buoyancy  forces;  or  forced  convection,  accomplished 
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Start 


Read: 

a,  k,  h,  radius,  number  of  points  in  fruit, 

ini- 

tial  product  temperature,  At,  total  cooling 

time 

print  frequency,  test  number,  specific  heat 

of 

air,  specific  heat  of  product,  total  weight 

of 

1 

product,  air  mass  flow  rate,  number  of  layers  in 

! 

bed,  bed  depth,  entering  air  temperature. 

Comoute  combinations  of  input  data  which  remain  constanl 


1 


=~l 


Compute  the  temperature  at  each  point  within  the  indi 
vidual  product  for  the  specified  time  period  and  the 
specified  layer  in  the  bed. 


Compute  the  mass-average  temperature  of  the  product  for 
the  specified  laver  in  the  bed. 


I  Compute  the  temperature  of  the  air  entering  the  next 
layer  of  the  bed  based  upon  the  amount  of  heat  lost  by 
the  product  in  the  preceding   layer. 


No 


Increase  time 
-,  by  one  time 
increment 


Print  time  and 
| all  temperatures 
1  for  that  time 
! oeriod. 


Figure  3.   Flow  diagram  for  the  digital  computer  program 
to  solve  deep  bed  heat  transfer  equations 
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RADIUS  RATIO 

Figure  9.   Magnitude  and  location  of  mass-average  temp- 
erature of  Marsh  grapefruit  cooled  in  agitated 
ice  water 
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Figure  10.   Unaccomplished  temperature  change  and  location 
of  mass-average  temperature  during  forced-air 
cooling 
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by  some  external  means  such  as  a  pump  or  fan.   The  principal 
resistance  to  heat  transfer  is  found  in  a  relatively 
stagnant  laminar  layer  adjacent  to  the  solid  surface.   In 
laminar  flow,  fluid  particles  move  only  in  a  direction 
parallel  to  the  solid  surface,  meaning  that  heat  transfer 
through  this  layer  must  be  accomplished  by  conduction  only. 
Adjacent  to  this  laminar  layer  is  usually  found  a  turbulent 
layer  of  slow  moving  particles  which  also  gives  resistance 
to  heat  transfer  between  the  main  bulk  fluid  and  the  solid 
surface.   These  two  regions  constitute  the  boundary  layer 
and  present  a  very  complex  problem  when  analyzed  mathematically. 
Thus,  the  study  of  convection  heat  transfer,  in  general,  has 
been  empirical  and  semi-empirical. 

The  rate  of  convection  heat  transfer  is  proportional  to 
the  difference  in  temperature  between  the  surface  and  the  main 
bulk  of  fluid  and  to  the  surface  area,  thus: 


q   =  h  A  (T   -  T  )  [65] 

■'C  S      oo 


The  proportionality  constant  (h)  is  called  the  convective 
heat  transfer  coefficient  and  is  determined  by  the  properties 
of  the  fluid,  the  nature  of  the  surface,  and  the  manner  and 
velocity  of  fluid  flow  past  the  surface. 

One  method  which  has  been  used  successfully  to  establish 
a  relationship  for  h  as  a  function  of  the  important  proper- 
ties of  the  system  is  dimensional  analysis.   This  method  of 
analysis  identifies  the  necessary  dimensionless  groupings  as 
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the  Nusselt  number  (hD/k) ,  the  Reynolds  number  (pDV/y) , 

the  Prandtl  nuniber  (cy/k)  ,  and  sometimes  a  shape  factor  (1/D) 

Numerous  experimental  and  semi-empirical  studies  have 
been  made  for  flow  over  various  shaped  objects,  utilizing 
these  dimensionless  parameters.   McAdams  (35)  recommended 
equation  [6  6]   for  calculating  the  average  convective  heat 
transfer  coefficient,  h,  for  spheres  heated  or  cooled  by  a 
gas  for  the  Reynolds  number  range  25  to  100,000. 


13s>  =  0.37   (VjzJDjO.e  [66] 

K  u 


where : 

D  =  diameter  of  sphere 

k  =  thermal  conductivity  of  the  gas 

V  =  velocity  of  fluid  stream 

p  =  density  of  fluid 

p  =  viscosity  of  fluid 

Several  equations  for  calculating  h  in  a  packed  bed 
have  been  recommended,  however,  they  give  quite  varied 
results.   Eckert  (16)  proposed  a  general  equation  for  h 
for  a  packed  bed  of  solid  particles: 


M=  0.8  (Re.)0'7  (^)1/3  [67] 
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where 

d  =  equivalent  spherical  diameter 

oVd 

Re .    =   Reynolds   number    (J )    where   V   is    the    mtersti- 

1  U 

tial  velocity  through  the  bed 
c  =  specific  heat  of  gas 

Daizo  and  Motoyuki  (13)  gave  the  heat  transfer  correla- 
tion as  : 


hD  =  1>2  [l=Jlzi±5/3]l/3  [2g6  ]V3     [68] 


where 

w  =  2-3r  +  5r5-2r6 

r-  (1-,)1/3 

£  =  voxd  fraction  of  packed  bed 

2 

G  =  superficial  mass  velocity  of  gas,  lb/hr-ft 

Lebedev  and  Petror-Denison  (32)  presented  a  correlation 

based  on  hydrodynamics : 


*£=    0.124    R   °'8    ?r0'33  [693 

k  e 


where 

Pr  =   Hg. 
k 

A  large  amount  of  experimental  information  on  heat 
transfer  in  packed  beds  has  been  analyzed  by  chemical  engi- 
neers (51)  to  arrive  at  the  following  empirical  correlations 
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^Pr2/S   =    0.91   Rem-°'51    *,    Rem<50  [70] 


JPr2/3=    0.61   Rem-°-41    *,      ReIft     >50  [71] 


where 


Re™  ~  a*  Mi,  =  modified  Reynolds  number 


_G_ 

~m  ""  A*yy 

A*  =  solid  particle  surface  area  per  unit  bed  volume 


The  quantity  \p  is  an  empirical  coefficient  that  depends  on 
the  particle  shape.  A  few  sample  values  of  ty  are  given  in 
Table  3(51). 

Table  3.   Particle  shape  Factors  for  packed-Bed  Correlation, 

Particle  shape  . 

Spheres  1.0  0 

Cylinders  0.91 

Flakes  0.86 

Raschig  rings  0.79 

Partition  rings  0.67 

Berl  saddles  0.80 

Barker  (6 )  has  made  an  extensive  review  of  literature 
on  heat  transfer  in  packed  and  fluidized  beds  and  has  com- 
piled a  bibliography  of  nearly  500  publications. 
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Smith  et  al.  (47)  describes  a  method  of  experimentally 
determining  the  convective  heat  transfer  coefficient  (h) 
and  the  thermal  diffusivity  (a)  for  regular  shaped  objects 
as  well  as  for  anomalous  shapes.   His  method  is  based  on  the 
solution  of  the  Fourier  conduction  equation  for  each  speci- 
fied geometry.   The  general  form  of  the  prediction  equation, 
for  boundary  conditions  of  finite  internal  and  surface  thermal 
resistance  is : 


0  =  f±     (Fo,  Bi,  Ld,  A,  3)  [72] 


where 


6  =  dimensionless  temperature  ratio 
Fo  =  — y  ,  Fourier  number 

Bi  =  — r-,  Biot  number 

L,  =  Distance  ratio  indicating  fractional  distance  from 

the  center 
A  =  a/£,  distance  ratio  used  to  define  the  geometry  of 

the  ellipsoidal  model  for  the  shape 
B  =  b/£ 
I   =   cnaracteristic  length  or  minimum  distance  within 

the  object  with  maximum  temperature  difference  during 

transient  period 
a  =  semimajor  axis  of  smaller  model  ellipse 
b  =  semimajor  axis  of  larger  model  ellipse 
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The  solution  of  the  heat  conduction  equation  for  a 
sphere  was  presented  in  the  previous  chapter  as: 

~>   sin  M  r 
co    sin  M   a  -  M   a  cos  M   a    .,2   n    ,   , 

0=2Z    (s Z—^-rr r^— )e~^nt   M   r     [21] 

..Ma-  sin  M  a  cos  M   a  n 

n=l   n  n        n 

M  is  the  separation  constant  arising  from  the  solution  of 
the  differential  equation  and  denotes  the  roots  of  the 
transcendental  equation. 


M  a  cos  M   a  =  U-Bi)  sin  14  a         [73] 
n        n  n 


Smith  et  ai.  (47)  suggested  that  equations  [21]  and  [73] 
for  various  geometries  could  be  represented  by 


8-  f2(A,B,Bi,Ls)e-Fo[f3(A<B'Bi>i2 


with  the  following  restrictions . 

Temperatures  in  the  solid  are  restricted  to  those  along 

the  characteristic  length;  the  geometry  may  be  described  by 

the  parameters  A  and  B;  only  the  first  term  of  the  infinite 

series  is  used;  and  M  is  used  as  a  parameter  for  any 

n  c 

specified  geometry.   For  a  specified  geometry,  equation  [74] 
could  be  written  as: 


9   =   Ce-M12F°  [75] 
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Then  a  geometry  index  G  is  defined  as  the  number  that, 

2 
when  multiplied  by  it  ,  gives  the  appropriate  value  for 

2 
the  coefficient  M..   in  equation  [75]  ,  for  the  geometry 

concerned  and  the  conditions  of  negligible  surface  thermal 

resistance. 

Experimental  values  for  G  ranged  from  0.25  to  1.0  0 

for  the  infinite  slab  and  the  sphere,  respectively.   The 

index  for  anamalons  shapes  is  described  by  equation  [76]. 


G  =  1/4  +  3/8A2  +  3/8B2  [76] 


For  conditions  of  negligible  surface  thermal  resis- 
tance (large  h)  and  considering  successive  values  of  0  and 
Fo,  Smith  et  al.  arrived  at  the  following  equation  for 
determining  thermal  diffusivity  (a) : 


2     In0,-ln09   . 
a  =  I—      [-—/_.   "  ]  177] 

M  2      *-!  "2 


For  negligible  surface  thermal  resistance: 


M12  =  Gu2  [7  8] 


For  conditions  other  than  negligible  surface  thermal 

resistance  (finite  h)  the  functional  relationship  between 

2 

M,   and  G  may  be  expressed  in  the  form: 

M  2  =  f3  (G,Bi)  [  79] 


56 

A  nomograph  for  evaluating  elements  of  this  equation  is 
given  in  Figure  11. 

After  a  has  been  determined  from  time-temperature  data 
taken  from  a  test  where  h  is  large,  such  as  an  agitated 
water  bath,  the  convection  heat  transfer  coefficient  (h) 
can  be  determined  from  the  following  equation  and  Figure  11. 


0    ,  2    In6„-ln80 


This  equation  is  analogous  to  equation  [77],  but  for  con- 

2 

ditions  where  h  is  finite  and  M,   is  no  longer  described 

2 
by  M1   =  GTT  . 

One  must  be  cautious  when  using  equations  [77]  and 
[80]  to  insure  that  the  boundary  and  initial  conditions  of 
the  original  conduction  equation  are  satisfied,  i.e.  uniform 
initial  temperature  distribution  and  constant  fluid  tempera- 
ture surrounding  the  object. 

Minh  et  al.  (36)  evidently  disregarded  the  boundary 
condition  of  constant  fluid  temperature  when  he  used 
equation  [30]  to  determine  h  at  three  depths  within  a  bed 
of  potatoes.   It  is  obvious  that  this  boundary  condition 
is  satisfied  only  for  the  layer  of  product  at  the  entrance 
of  the  bed.   Their  values  for  h  decreased  as  distance  into 
bed  increased,  with  the  entrance  value  being  approximately 
three  times  the  exit  value.   They  developed  the  following 
equation  for  white  potatoes  at  the  entrance  to  the  bed: 


■si 

<i  ■ 

7  | 

(. 

5: 
A  - 


57 


M    Slab 


1.0  r 
0'  - 
0.8 1 
07 

06 

05 

01 


0  2 


■ 


0.10 

oost- 

0.08?- 

0  0:'- 
0  00- 
0.D     - 

c.t-.,- 


o.o:;- 


-2.0 


-20 


-30 

-<0 

-50 

-60 
-70 
-80 
-30 
(-100 


-025 


hO  30 


r-0.40 


0.50  £ 
E 


Lf.  Cjl.- 


:-0.60 


-0.70 


-0.80 
I 
^0.90. 


-/....  -f- 


1.00 


oo;r- 
t 


O.Olf- 


Figure  11.   Nomograph  for  evaluating  elements  of  the 
equation:  Hi 


=  f3(G,3i) 


58 


*A=  0.564  (Re)0"538  (Pr)1/3        [8l] 
is. 


Bennett  et  al.  (9)  conducted  tests  on  the  cooling  of 

citrus  fruits  in  bulk  and  developed  a  prediction  equation 

for  what  they  called  an  effective  h  (h  ,.,.)  as  indicated  in 
J  efr 

Newton's  equation: 


h  ^  =  Q/(AAt)  [82] 


where  for  this  equation,  Q  is  the  total  heat  removed  from 
the  mass  of  fruit  per  hour,  A  is  the  total  surface  area  of 
all  the  fruit  based  on  average  diameter,  and  At  is  the 
average  difference  between  temperature  on  the  surface  of  the 
fruit  and  that  of  the  air.   They  did  not  attempt  to  account 
for  any  variation  of  h  as  a  function  of  position  within  the 
bed.   By  averaging  all  the  values  obtained  at  various  points 
in  the  bed,  they  obtained  the  following  heat  transfer  corre- 
lations : 


^=C(Re.)n  [83] 


where  C  and  n  arc  constants  that  were  experimentally  evaluated 
as  indicated  in  Table  4.   The  correlation  coefficient  is 
denoted  by  r.   Their  reported  correlation  coefficients  are 
indicative  of  the  difficultv  encountered  in  obtaining 
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reliable  data.   In  the  tests  from  which  these  data  were 
obtained  the  air  entering  the  bed  was  not  held  constant, 
thus  making  it  difficult  to  analyze  the  results  and  prevents 
the  determination  of  h  by  the  method  described  by  Smith 
ct  al.  (46) . 


Table  4.   Coefficients  and  Exponents  for  Emperical  Heat 
Transfer  Correlations 


C  N  r 

Orange                   6.95  0.297  0.610 

Grapefruit              40.9  0.153  0.204 

Tangelo                  6.30  0.314  0.464 

Combined                 5.05  0.333  0.536 


Development  of  Method 

A  method  similar  to  that  described  by  Bennett  et  al.  (9  ) 
utilizing  Newton's  equation  was  used  in  an  attempt  to 
investigate  the  variation  of  h  with  respect  to  depth  in  a 
bed  of  oranges.   However,  preliminary  results  indicated  the 
precision  of  the  measurements  was  not  sufficient  to  detect 
any  variation  in  h  using  this  method.   It  is  very  difficult 
to  obtain  high  level  precision  using  this  method  since  the 
placement  of  the  thermocouples  is  very  critical  in  addition 
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to  the  normal  variation  which  occurs  in  a  bed  of  irregularly- 
shaped  particles.   The  effect  of  this  irregularity  on  the 
determination  of  h  is  greater  for  this  method  than  it  is  for 
the  method  described  by  Smith  et  al.   Bennett  et  al . 
utilized  surface  temperatures  which  are  more  sensitive  to 
local  variations  in  the  surrounding  fluid  and  are  in  general 
more  difficult  to  measure  than  internal  temperatures.   One 
notable   advantage  of  the  method  described  by  Smith  et  al.  is 
that,  theoretically,  it  is  independent  of  the  location  of 
the  thermocouple  and  thus  eliminates  that  source  of  error. 
The  method  for  determining  h,  described  by  Smith  et  al .  and 
previously  outlined  in  this  section, was  used  to  evaluate  h 
for  the  entrance  layer  of  the  bed.   The  values  for  h  were 
determined  for  beds  of  oranges  and  grapefruit  as  a  function 
of  superficial  air  velocity  (V)  entering  the  bed.   The 
results  are  presented  in  Figures  20  and  21  (Data  and  Results 
section)  in  terms  of  a  Nusselt  number  versus  Reynolds  number 
plot.   Since  these  h  values  were  determined  only  for  the 
bottom  layer  of  fruit,  it  is  expected  that  they  do  not 
represent  the  best  average  value  for  the  entire  bed. 

To  overcome  this  deficiency,  another  method  for  eval- 
uating h  was  investigated.   Using  the  h  values  determined 
for  the  bottom  layer,  the  computer  model  for  deep  bed  cool- 
ing, as  described  by  Figure  8,  was  run  for  each  air  flow 
rate.   Mass  average  temperatures  from  the  model  were  com- 
pared with  those  from  the  experimental  data  for  each  test. 
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From  these  comparisons,  h  values  were  estimated  so  that 
at  least  one  value  would  be  lower  and  one  value  would  be 
higher  than  the  value  which  gave  the  best  fit  for  the 
experimental  data.   Then  these  values  were  used  to  generate 
more  data  from  the  computer  model.   A  comparison  of  the 
experimental  and  theoretical  mass  average  temperature 
of  the  entire  bed,  determined  at  approximately  "one  half 
cooling  time,"  was  used  as  the  criterion  for  determining 
the  h  value  which  resulted  in  the  best  fit  to  the  experi- 
mental data.   A  one  half  cooling  time  is  the  time  required 
to  reduce  the  fruit  temperature  by  one  half  of  the  difference 
between  the  initial  product  temperature  and  the  temperature 
of  the  cooling  air.   These  results  are  also  shown  in  Figures 
20,  21,  and  22  (Data  and  Results  section). 

A  linear  regression,  by  least  squares,  of  the  log  N 
versus  log  Re  was  used  to  evaluate  the  constants  c  and  n 
in  equation  [84 ] . 


Nu  =  c(Re)n  lo-.j 


where : 

Nu  =  -— ,  Nusselt  number 

Re  =  — —  ,  Reynolds  number 

A  separate  regression  was  made  for  each  group  of  data  as 
indicated  in  Figures  20,  21,  and  22  and  Table  5.   Theoretically, 
oranges  and  grapefruit  should  be  governed  by  the  same  equa- 
tion.  However,  some  variation  would  be  expected  since  a 
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homogeneous  sphere  is  only  an  approximation  of  a  grapefruit 
and  an  oranqe. 


Table  5.   Coefficients,  Exponents  and  Linear  Correlation 
Coefficients  for  Equation  [04]  Determined  for 
Forced-Air  Cooling  in  Bulk  Containers. 


c  n  r 

Orange 
(Bottom  layer)  0.859  0.521       0.964 

Orange 
(Average  for  bed)         0.374  0.620       0.991 

Grapefruit 
(Bottom  layer)  7.17  0.272       0.832 

Grapefruit 
(Average  for  bed)         1.04  0.482       0.987 

Oranges  and  Grapefruit 
Combined 
(Average  for  bed)         0.607  0.555       0.977 


Experimental  Evaluation  of  Thermal  Diffusivity 

The  thermal  diffusivity  (a)  was  determined  experimentally 
for  oranges  and  grapefruit  by  the  method  described  by  Smith 
et  al.   (47).   The  fruit  were  cooled  in  an  agitated  water  bath 
with  a  thermocouple  placed  at  the  approximate  mass  average 
point.   Theoretically,  the  results  are  independent  of  the 
thermocouple  position.   The  temperatures  were  recorded  by 
a  six-point  potentiometer  with  a  print  frequency  of  three 
seconds.   Equation  [77]  was  used  to  calculate  a,  with  the 
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loge8l"  loge62 
term z determined  from  a  semi-log  plot  of 

rl  r2 

the  experimental  data.   M.   was  calculated  using  equation 

[76]  and  [78].   Oranges  were  considered  perfect  spheres, 

while  G  for  grapefruit  was  calculated  from  equation  [77], 

based  on  the  equatorial  and  polar  diameters  of  each  fruit 

tested.   The  thermal  diffusivity  (a)  for  oranges  was  found 

to  be  0.00437  ft  /hr  with  a  standard  deviation  of  0.00010 

based  on  five  tests,  while  a  for  grapefruit  was  found  to 

be  0.00362  ft2/hr  with  a  standard  deviation  of  0.00015 

based  on  five  tests. 


Design  and  Description  of  Experimental  Facilities 

A  number  of  investigators  have  conducted  research  on 
precooling  using  forced-air  cooling  equipment.   Most  of  these 
investigators  emphasized  the  response  of  the  product  to  cool- 
ing.  A  lesser  number  were  concerned  with  the  variable  factors 
which  influence  the  design  and  operation  of  cooling  equipment. 

A  forced-air  cooler  designed  specifically  for  research, 
with  a  capability  of  controlling  all  of  the  important  vari- 
ables relating  to  cooling  biological  materials,  was  con- 
structed for  this  study.   Air  velocity,  air  temperature, 
relative  humidity,  and  product  arrangement  were  the  variables 
that  could  be  controlled.   Temperature  distribution  within 
the  individual  product  and  within  the  bulk  load,  static 
pressure  loss  through  the  product,  and  product  moisture 
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loss  were  among  the  parameters  that  could  be  measured 
during  a  test. 

Basic  Construction  of  Precooler 

Figure  12  shows  the  general  location  of  the  components 
and  the  air  circulation  diagram  for  the  cooling  facility. 
The  overall  length  of  the  cooler  is  25  feet  and  the  outside 
dimensions  of  the  product  chamber  section  are  8x8x12  feet. 
The  product  chamber  was  designed  to  accommodate  a  product 
load  which  could  have  dimensions  of  up  to  4x4  feet  in  cross- 
section  and  8  feet  high.   This  would  allow  experimental 
units  of  pallet  boxes  or  pallet  loads.   The  basic  wall  con- 
struction is  2x4-inch  studding  with  1/2-inch  plywood  on  each 
side.   Insulation  material  consisted  of  rigid  polyurethane 
foamed  in  place  between  the  plywood  sheets.   The  outside 
of  the  cooler  was  sealed  with  caulking  and  epoxy  paint. 
The  complete  structure  is  built  on  skids  and  can  be  moved 
when  necessary. 

Refrigeration  Components  and  Controls 

Figure  13  shows  each  component  of  the  refrigeration 
system.   The  condensing  unit  is  water  cooled  and  is  preferred 
over  an  air-cooled  unit  when  precise  control  is  necessary, 
since  changes  in  ambient  conditions  have  little  effect  on 
its  performance.   The  two-cylinder  compressor  is  belt  driven 
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by  a  7-1/2-hp  motor  and  can  accommodate  up  to  20-hp  if 
more  capacity  is  needed.   The  compressor  is  equipped  with 
an  automatic  cylinder  unloader  which  inactivates  one  of 
the  cylinders  when  the  suction  pressure  drops  below  an 
adjustable  preset  value.   A  belt-driven  compressor  is  de- 
sirable when  the  system  is  to  be  used  for  various  applica- 
tions where  the  evaporator  temperatures  may  vary  over  a 
wide  range.   The  following  example  will  illustrate  this 
point. 

Example  illustrating  capacity  of  syster 

Using  Refrigerant  502  with  a  compressor  rpm  of  520, 
evaporator  temperature  of  30  F  and  a  condenser  temperature 
of  100  F,  7  tons  of  cooling  will  be  produced,  requiring 
7-1/2-hp.   If  the  evaporator  temperature  is  increased  to 
55  F,  with  other  conditions  unchanged,  approximately  11 
tons  of  cooling  will  be  produced  requiring  6-1/2-hp.   When 
the  evaporator  temperature  is  decreased  to  -20  F,  approxi- 
mately 2-1/2  tons  of  cooling  will  be  produced,  requiring 
6  hp .   Since  cooling  capacity  and  horsepower  vary  with 
compressor  rpm,  it  may  be  necessary  to  change  the  compressor 
rpm  in  order  to  obtain  the  desired  cooling  capacity  at  a 
given  evaporator  temperature  or  to  adjust  the  horsepower 
required  at  a  certain  evaporator  temperature. 


Components  used  for  control 


The  refrigeration  system  is  designed  to  allow  continuous 
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operation  of  the  compressor.   Intermittent  operation  would 
result  in  poor  temperature  control.   Air  temperature  control 
is  accomplished  by  adjusting  the  evaporator  pressure  with  a 
motorized  pressure  regulating  valve  on  the  suction  side  of 
the  evaporator.   This  valve  is  sometimes  called  a  back- 
pressure regulator  or  an  evaporator  pressure  regulator. 

Since  sometimes  the  cooling  load  varies  over  a  wide 
range  during  a  single  test  while  the  compressor  is  operat- 
ing continuously,  it  is  necessary  to  have  a  hot-gas  bypass 
valve  in  addition  to  the  cylinder  unloader  to  provide  for 
additional  capacity  reduction.   As  indicated  in  Figure  13, 
this  valve  recirculates  hot  gas  from  the  compressor  back 
through  the  evaporator  when  the  suction  pressure  drops  below 
an  adjustable,  preset  level.   The  motorized  pressure 
regulating  valve  is  operated  by  an  automatic  controller 
which  will  be  discussed  in  the  next  section. 

Evaporator  coils 

The  direct  expansion  evaporator  coils  shown  in  Figure  13 
are  standard  models  without  housing  or  fan.   They  have  wide- 
spaced  fins  (1/4  inch)  so  that  a  thin  layer  of  ice  collected 
during  a  test  will  not  significantly  affect  the  air  flow 
rate.   Both  evaporators  are  equipped  with  defrosting  heaters 
which  can  be  actuated  manually.   The  heaters  can  also  be 
used  to  provide  additional  heat  load  on  the  refrigeration 
system  which  is  sometimes  desirable  in  order  to  achieve  better 
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control.   Refrigerant  502  is  used  in  this  system  since  it 
allows  a  wide  range  of  operating  temperatures  without 
changing  the  thermal  expansion  valves  and  other  control 
components . 

Air  Temperature  and  Humidity 

The  automatic  control  of  air  temperature  entering 
the  product  chamber  is  accomplished  by  the  use  of  a  propor- 
tional controller  with  reset  and  rate  action,  which  regulates 
the  motorized  valve  on  the  evaporator.   The  primary  sensing 
element  for  the  automatic  controller  is  a  thermocouple 
which  is  placed  in  the  air  stream  at  the  entrance  of  the 
bed.   The  manipulated  variable  is  the  evaporator  pressure 
which  is  regulated  by  the  motorized  valve. 

For  tests  during  which  the  inlet  air  temperature  is  to 
remain  constant,  the  cooling  load  at  the  beginning  of  the 
test  will  be  much  greater  than  the  average  cooling  load. 
Thus,  it  is  desirable  for  a  research  facility  to  have  a 
refrigeration  system  large  enough  to  handle  a  temporary 
high  cooling  load  for  products  that  cool  extremely  fast. 

Since  the  air  is  recirculated,  the  dew  point  of  the 
cooling  air  will  approach  the  evaporator  temperature  unless 
moisture  is  added.   In  some  cases  an  undesirable  drying 
effect  is  caused  by  a  large  water  vapor  pressure  difference 
between  air  and  the  product  surface.   In  order  to  minimize 
this  effect,  the  cooler  is  equipped  with  a  steam  humidifier 
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which  injects  steam  into  the  air  after  it  passes  through 
the  evaporator  coils.   Adding  steam  to  the  air  in  this 
manner  results  primarily  in  a  latent  heat  increase,  and 
very  little  sensible  heating  occurs.   For  example,  when 
air  at  30  F  and  70%  RH  is  saturated  by  adding  212  F  steam, 
the  air  temperature  is  increased  less  than  0.4  F.   This 
procedure  presents  no  problem  in  controlling  the  air 
temperature  since  the  temperature  sensor  is  downstream 
from  the  humidifier.   Using  large  evaporators,  thus  keep- 
ing the  temperature  difference  between  the  coil  and  the 
air  small,  will  result  in  higher  humidities.   However, 
the  humidity  will  always  be  below  100%  if  moisture  is  not 
added  back  to  the  air  after  passing  through  the  coils.   Even 
when  the  air  entering  the  product  chamber  is  saturated, 
there  is  still  a  potential  for  drying  in  a  product  bed 
because  the  air  is  heated  by  the  preceding  layer  of  product. 

In  order  to  determine  the  amount  of  moisture  loss  during 
cooling,  a  load  cell  type  of  scale  was  mounted  so  that  the 
entire  product  bin  could  be  suspended  and  weighed  before 
and  after  each  test. 

Air  Flow 

The  location  of  the  air  handling  unit  and  air  flow 
patterns  are  shown  in  Figure  12.  The  air  is  circulated 
with  a  backward  curved  centrifugal  blower  having  a  capacity 
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of  16,000  cfm  when  operating  against  3  inches  H_0  static 
pressure.   The  backward  curved  vanes  produce  a  relatively 
high  operating  pressure  and  will  provide  approach  velocities 
of  up  to  1000  fpm  entering  the  product  chamber. 

The  blower  is  driven  by  a  hydraulic  motor  controlled 
by  a  variable  speed  hydraulic  transmission  which  allows 
the  air  flow  rate  to  be  regulated  by  the  rpm  of  the  blower 
rather  than  by  dampers.   The  transmission  and  electric 
motor  are  located  outside  of  the  cooler,  which  minimizes 
the  heat  produced  inside  the  cooler.   The  blower  was 
placed  downstream  from  the  evaporator  coils  so  that  the 
temperature  of  the  air  through  the  blower  could  remain 
constant  throughout  a  test.   If  the  blower  were  placed 
upstream  from  the  evaporators,  the  temperature  of  air 
moving  through  the  fan  would  vary  during  the  test  causing 
a  variation  in  the  velocity  through  the  product  chamber. 

The  air  enters  a  large  plenum  and  passes  through  an 
air-straightening  section  beneath  the  product  chamber  to 
ensure  that  the  air  entering  the  product  is  uniformly 
distributed  as  it  enters  the  product  chamber.   In  conjunc- 
tion with  the  air-straightening  vanes,  small  patches  of 
screen  wire  were  added  to  produce  resistance  in  areas  which 
have  a  high  velocity.   The  velocity  profile  of  the  air 
entering  the  product  chamber  did  not  vary  more  than  5%  when 
measured  with  a  hot-wire  anemometer  probe  at  6-inch  intervals 
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The  air  flow  rate  is  measured  to  within  -  2  % 
through  the  use  of  an  "annubar"  flow  element  installed 
in  a  1-foot-diameter  PVC  pipe.   The  flow  rate  is  determined 
as  a  function  of  the  difference  between  velocity  and  static 
pressure.   The  pressure  difference  is  measured  with  an 
electronic  differential  pressure  indicator  which  has  a  reso- 
lution of  0.001  mm  Hg. 

Reheat  Section 

The  reheat  section,  as  shown  in  Figure  12,  is  a  self- 
contained  unit  which  may  be  detached  from  the  main  struc- 
ture or  placed  adjacent  to  the  product  chamber  so  that  a 
container  load  of  product  can  be  rolled  out  and  reheated 
for  repeated  tests  on  the  same  product  load.   It  is  also 
used  to  heat  the  product  to  a  predetermined  uniform  temper- 
ature before  the  initial  test.   It  is  equipped  with  a 
blower,  heater,  temperature  controller,  and  humidifier. 
Since  the  product  is  heated  outside  of  the  cooling  chamber, 
the  cooler  can  operate  continuously  and  maintain  the 
desired  temperature  while  the  product  is  being  heated. 
The  primary  source  of  heat  for  the  heating  section  is  a  steam- 
heated  coil  with  secondary  electric  heaters  used  for  fine 
control.   The  electric  heaters  are  regulated  by  a  proportional 
automatic  controller. 

Experimental  Test  Procedures 
Ten  tests  were  conducted  on  a  340-pound  bulk  load  of 
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No.  100  Valencia  oranges  with  superficial  air  velocities 
ranging  from  25  to  40  0  fpm.   These  tests  were  conducted 
on  two  different  bulk  loads,  five  tests  each,  in  order  to 
get  some  information  on  the  variability  among  bulk,  loads 
of  oranges.   Five  tests  were  conducted  on  a  274-pound  bulk 
load  of  No.  23  Marsh  grapefruit  with  superficial  air 
velocities  ranging  from  15  to  4  00  fpm. 

Each  cooling  test  was  conducted  using  eight  cartons 
(4/5  bushel  each)  of  sized  fruit  which  were  weighed, 
counted,  and  placed  into  the  bulk  container.   This  con- 
tainer was  2x2x4  feet  deep  with  an  expanded  metal  floor. 
Thermocouples  were  placed  in  eight  fruit  which  were 
positioned,  two  each,  at  four  levels  as  indicated  in 
Figure  14 .   These  thermocouples  were  placed  at  the  mass 
average  temperature  point  and  at  the  surface.   Center 
temperatures  were  not  measured  since  citrus  fruit  usually 
have  a  void  space  at  the  center.   The  thermocouples  were 
constructed  from  36-gage,  insulated  copper-constantan  wire. 

Figure  15  depicts  the  method  used  in  placing  the  thermo- 
couples in  the  product.   The  insulated  thermocouple  wire 
was  first  inserted  into  a  hypodermic  needle  and  then  injected 
into  the  product  to  the  desired  location.   Then  the  needle 
was  carefully  withdrawn  leaving  the  thermocouple  in  place. 
The  thermocouple  was  then  sealed  at  the  surface  of  the 
product.   A  reference  mark,  which  always  remains  outside 
of  the  product,  was  placed  on  all  thermocouple  wires  to 
indicate  the  thermocouple  position  in  case  of  slippage. 
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Figure  14.   Bulk  load  of  fruit  showing  the  four  levels  at 

which  product  temperature  was  measured.   Dimen- 
sions A,  B,  C,  and  D  are  indicated  for  each  test 
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Bennett  et  al.  (8)  have  suggested  that  errors  due  to  con- 
ductance along  the  thermocouple  wire  from  the  surrounding 
cold  air  could  be  minimized  if  the  thermocouple  was  always 
placed  at  least  halfway  through  the  product  in  order  to 
produce  a  reverse  temperature  gradient  in  the  thermocouple 
wire.   This  method  was  used  for  temperature  measurement  at 
the  mass  average  point,  but  preliminary  tests  on  plastic 
spheres  indicated  that  a  considerable  error  could  result 
when  surface  temperatures  were  measured  by  inserting  the 
thermocouple  through  the  sphere  and  outward  to  the  surface. 
The  surface  temperatures  measured  in  this  manner  were  always 
higher  than  those  calculated  by  using  the  center  or  mass 
average  temperatures  and  the  theoretical  solution  for 
conduction  in  a  sphere.   This  is  probably  due  to  the  fact 
that  the  thermocouple  is  located  in  an  outer  layer  of  the 
sphere  and  not  on  the  surface.   The  large  temperature  gradients 
in  the  outer  layers  of  the  sphere  also  could  cause  a  con- 
duction error  along  the  wire  near  the  junction.   In  an 
effort  to  improve  this  measurement,  the  conductors  were 
twisted  together  for  approximately  1/4  inch  and  imbedded 
along  the  surface  of  the  fruit. 

All  thermocouple  junctions  were  coated  with  a  thin  film 
of  plastic  to  prevent  possible  corrosive  or  electrolytic 
action.   The  calibration  of  all  thermocouples  was   checked 
at  32  F  and  80  F  and  was   found  to  have  an  error  of  less 
than  0.5  F. 
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Temperature  was  recorded  on  a  multi-point  potentiom- 
eter which  read  each  point  at  two-minute  intervals  through- 
out the  test.   Each  product  load  was  brought  to  a  uniform 
temperature  just  prior  to  placing  the  container  into  the 
cooling  chamber. 

Data  and  Results 
Numerical  Solution  of  Deep  Bed  Cooling 

Figures  16  through  19  show  typical  output  data  from  the 
numerical  solution  in  graphical  form.   The  computer  program 
and  sample  output  data  are  shown  in  Appendix  B.   It  should 
be  noted  that  only  a  small  portion  of  the  output  can  be 
indicated  on  any  one  figure  and  have  the  graphs  remain  legible, 
It  is  possible  to  produce  a  graph  of  the  temperature  at  any 
point  within  the  individual  fruit  and  at  any  depth  within 
the  bed,  as  a  function  of  time.   The  computer  program 
provides  the  capability  to  vary  the  important  thermal  and 
physical  properties  of  the  cooling  system  such  as:   the  size 
of  product,  thermal  properties  of  product,  air  velocity, 
air  temperature,  and  dimensions  of  the  bed.   A  complete  list 
of  these  variables  is  shown  in  Appendix  B  as  input  data. 

The  graphical  results  shown  in  this  section  are  only 
samples  from  the  complete  range  of  tests.   Additional 
results  from  the  numerical  model  are  shown  in  Appendix  C 
in  the  form  of  a  comparison  between  the  numerical  model  and 
experimental  data. 
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Figure  16.   Theoretical  mass-average  temperature  response  of 
oranges  to  forced-air  cooling  as  predicted  by- 
numerical  model 
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Figure  17.   Theoretical  mass-average  and  surface  temperature 
response  of  oranges  to  forced-air  cooling  as 
predicted  by  numerical  model 
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Theoretical  mass-average  and  surface  temperature 
response  of  oranges  to  forced-air  cooling  as 
predicted  by  numerical  model 
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Figure  19.   Theoretical  mass-average  and  surface  temperature 
response  of  grapefruit  to  forced-air  cooling  as 
predicted  by  numerical  model 
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Convective  Heat  Transfer  Coefficient  (h) 

A  relationship  for  h  in  a  deep  bed  of  citrus  fruit  was 
determined  in  terms  of  Nusselt  number  versus  Reynolds  num- 
ber as  indicated  by  the  following  equation. 


Nu  =  c(Re)n 


where 

Nu  =  hD/k,  Nusselt  number 
Re  =   VD/y,  Reynolds  number 

A  separate  regression  was  made  for  each  group  of  data  as 
indicated  in  Figures  20,  21,  22,  and  Table  5. 

Thermal  Diffusivity  (a) 

Thermal  diffusivity  was  determined  experimentally  for 
oranges  and  grapefruit  using  the  method  described  by  Smith 
et  al.  (47) ,  which  utilizes  time-temperature  data  from  fruit 
cooled  in  an  agitated  water  bath.   The  thermal  diffusivity  (a) 
for  oranges  was  found  to  be  0.00437  ft  /hr  with  a  standard 

deviation  of  0.00010,  based  on  five  tests,  while  a  for 

2 

grapefruit  was  found  to  be  0.00362  ft  /hr  with  a  standard 

deviation  of  0.00015,  based  on  five  tests. 

Experimental  Deep  Bed  Cooling  Tests 

The  experimental  temperature  response  from  forced~air 
cooling  tests  conducted  on  bulk  loads  of  Valencia  oranges 
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Figure  20.   Linear  regression  relating  Nusselt  number  and 
Reynolds  number  for  forced-air  cooling  of 
oranges  in  bulk  loads 
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Linear  regression  relating  Nusselt  number  and 
Reynolds  number  for  forced-air  cooling  of 
grapefruit  in  bulk  loads 
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Figure    22.      Linear   regression   relating   Nusselt   number   and 
Reynolds   number    for   forced-air   cooling   of 
oranges    and   grapefruit   in   bulk    loads 
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and  Marsh  grapefruit  are  shown  in  Figure  23  through  26  • 
Temperature  was  measured  at  the  surface  and  mass  average 
temperature  point  for  fruit  positioned  in  the  bulk  load  as 
indicated  in  Figure  14.   The  dimensions  A,  B,  C,  and  D 
indicated  in  Figure  14  are  specified  on  each  graph. 

The  results  presented  in  this  section  (Figures  23 
through  26)  do  not  include  all  of  the  experimental  tests. 
In  all,  there  were  fifteen  tests  run  on  oranges  and  grape- 
fruit with  superficial  air  velocities  ranging  from  15  to 
400  fpm.   The  complete  results  of  these  tests  are  contained 
in  Appendix  C,  in  the  form  of  a  comparison  between  the 
numerical  model  and  experimental  data. 

Table  6  gives  the  thermal  and  physical  properties  of 
the  Valencia  oranges  and  Marsh  grapefruit  used  for  the 
experimental  tests.   All  of  these  properties  were  deter- 
mined experimentally  or  calculated  using  experimental  data, 
except  specific  heat  (c) ,  which  was  taken  from  reference  (2) 

Discussion 


Comparison  of  Laplace  Transform  Solution  for  the  Simplified 
Model  and  Numerical  Solution  for  the  Realistic  Model 


In  order  to  obtain  an  analytical  solution  for  cooling  of 
a  deep  bed  of  citrus  fruit,  a  simplifying  assumption  of  no 
temperature  gradient  within  the  individual  fruit  was  made. 
In  the  development  of  the  numerical  solution  this  assumption 
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Figure  23.   Experimental  mass-average  and  surface  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  24.   Experimental  mass-average  temperature  response 
of  oranges  to  forced-air  cooling 
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Figure  25.   Experimental  mass-average  and  surface  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  26.   Experimental  mass-average  temperature  response 
of  grapefruit  to  forced-air  cooling 
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Table  6.   Thermal  and  Physical  Properties  of  Valencia  Oranges 
and  Marsh  Grapefruit  Used  in  Cooling  Tests. 


Total  weight,  lb 

Total  number  of  fruit 

3 


Valencia        Marsh 
oranges       grapefruit 


Bulk  density,  lb/ft 
Specific  gravity 
Void  space  fraction 

Equivalent  diameter,  ft 

2    3 
Surface  area  per  volume,  ft  /ft 

Thermal  conductivity,  BTU/hr-ft-°F   0.248 

Specific  heat,  BTU/lb-°F 

2 

Thermal  diffusivity,  ft  /hr 


340 

274 

800 

184 

37.8 

31.0 

1.016 

0.843 

0.406 

0.400 

0.234 

0.380 

15.3 

9.42 

0.248 

0.164 

0.90 

0.91 

0.00437 

0.00342 

was  lifted,  thus,  accounting  for  conduction  heat  transfer  with- 
in the  individual  fruit.   For  very  slow  cooling  rates, 
corresponding  to  low  air  flow  rates ,  both  solutions  should 
agree  since  there  will  be  a  negligible  temperature  gradient 
within  the  individual  fruit.   Figure  27  shows  that  for  an 
air  flow  rate  of  approximately  25  fpm  there  was  no  detectable 
difference  between  the  two  solutions. 
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Laplace  transform  solution 
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Figure  27.  Comparison  of  Laplace  transform  solution  of 
simplified  model  with  numerical  solution  of 
realistic  model  for  a  slow  cooling  rate  (25fpm) 
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This  comparison  indicates  that  the  simplified  model 
is  valid  for  low  air  flow  rates  and  also  provides  a  check 
on  the  numerical  solution.   A  similar  comparison  was  made 
in  the  preceding   chapter  (Figure  4)  which  provided  a  check 
on  the  numerical  solution  for  the  case  where  conduction 
within  the  fruit  was  important. 

In  order  to  demonstrate  the  effect  of  significant  tem- 
perature gradients  within  the  individual  fruit  on  the 
cooling  rate,  conditions  corresponding  to  a  faster  cooling 
rate  were  used.   Figure  28  gives  this  comparison  for  an  air 
velocity  of  approximately  400  fpm.   Thus,  it  is  seen  that 
the  results  are  considerably  different  when  the  effect  of 
temperature  gradients  within  the  individual  fruit  is 
considered. 

The  four  curves  on  Figures  27  and  28  indicate  the 
mass-average  temperature  at  four  levels  within  the  bed  of 
fruit  as  illustrated  in  Figure  14. 

The  accuracy  of  the  numerical  solution  is  obviously 
dependent  upon  the  number  of  increments  within  the  individ- 
ual fruit,  the  number  of  increments  within  the  bed  and  the 
number  of  time  increments.   As  explained  in  the  preceding 
chapter,  the  ratio  of  the  number  of  time  increments  to  the 
number  of  increments  within  the  fruit  has  a  lower  limit  in 
order  to  insure  stability  of  the  numerical  solution.   There- 
fore, the  number  of  increments  within  the  fruit  is  chosen 
first. 

A  series  of  tests  was  run  using  the  computer  program  to 
determine  the  effect  of  the  number  of  segments  in  the  fruit 
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Figure  28.  Comparison  of  Laplace  transform  solution  of 
simplified  model  with  numerical  solution  of 
realistic  model  for  a  fast  cooling  rate  (400fpm) 
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and  the  number  of  layers  in  the  bed,  on  the  fruit  tempera- 
ture.  The  effect  on  the  fruit  temperature  with  respect  to 
position  in  fruit,  position  in  bed,  and  time  was  studied 
using  the  computer  output.   The  largest  change  in  fruit 
temperature  as  a  function  of  number  of  layers  in  the  bed 
and  number  of  segments  in  the  fruit  was  determined  and  is 
presented  in  Table  7.   Using  this  information  and  consider- 
ing the  cost  of  computer  time,  the  following  number  of 
increments  was  chosen:   segments  in  fruit  -  10,  layers  in 
bed  -  5  per  foot,  and  time  increment  -  0.005  hours.   A 
decrease  in  the  time  increment  below  this  value  did  not 
indicate  a  detectable  increase  in  the  accuracy.   The 
results  shown  in  Figures  27  and  2  8  and  ail  other  results 
from  the  prediction  model  were  obtained  using  these  incre- 
ments .   Using  these  increments  for  a  cooling  test  on  a 
2-foot-deep  bed,  the  computer  time  was  less  than  2  seconds 
per  hour  of  cooling  time  when  run  on  an  IBM  370/165  model 
digital  computer. 

Convective  Heat  Transfer  Coefficient  (h) 

A  review  of  the  prediction  equations  for  h  in  packed 
beds  (pages  50-52)  indicates  very  little  consistency  among 
the  numerous  relationships.   The  problem  of  selecting  the 
correct  equation  has  perplexed  many  investigators  who  have 
studied  heat  transfer  in  packed  beds.   In  many  cases  it 


96 


Table  7.   The  Effect  of  Increment  Size  on  the  Accuracy  of 
the  Numerical  Solution. 


No.  of  Layers  Product  Surface  Temperature  for 

in  Bed  Top  Layer   (°F) 


5  77.4 

10  78.0 

15  78.2 


No.  of  Segments        Product  Mass  Average  Temperature 
in  Fruit  for  Bottom  Layer  (°F) 


5  66.7 

10  67.1 

15  67.3 


appears  that  too  much  effort  was  expended  on  one  facet  of 
the  problem  without  due  concern  for  how  the  results  would 
be  used  in  the  solution  of  the  total  problem. 

In  this  study  the  problem  of  selecting  the  correct 
equation  for  h  was  avoided  by  experimentally  determining 
the  relationship  for  h  utilizing  the  prediction  model 
developed  in  this  study.   This  relationship,  which  was 
presented  in  terms  of  the  Reynolds  number  and  Nusselt  num- 
ber (equation  [84])/  compares  favorably  with  relationships 
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used  by  other  investigators  of  heat  transfer  in  beds  of 
biological  materials.   The  experimentally  determined  values 
for  n  (equation  [84])  for  oranges  and  grapefruit  were  0.62 
and  0.48  respectively  as  indicated  in  Table  5.   Bakker- 
Arkema  (5)  and  Huang  et  al.   (29),  in  their  work  with  heat 
transfer  in  beds  of  beets  and  onions,  used  a  relationship 
with  an  n  equal  to  0.58  and  obtained  satisfactory  results. 

A  relationship  for  h,  valid  only  for  the  bottom  layer 
of  the  bed,  was  determined  using  a  completely  different 
method  as  described  by  Smith  et  al.   (27).   Figure  20  shows 
a  comparison  of  the  relationship  for  the  bottom  layer  and 
the  average  for  the  bed,  indicating  a  close  agreement. 

Theoretically,  h  for  oranges  and  grapefruit  would  be 
governed  by  the  same  equation  if  both  were  perfect  spheres. 
Thus,  a  linear  regression  was  made  for  all  experimental 
tests  conducted  on  oranges  and  grapefruit  combined.   Figure 
20  shows  that  a  good  correlation  was  obtained  with  a  corre- 
lation coefficient  (r)  equal  to  0.977. 

The  fact  that  h  determined  for  the  bottom  layer  only 
and  for  the  entire  bed  has  similar  values  indicates  that 
there  was  no  significant  variation  of  h  through  the  bed. 
The  value  of  h,  of  course,  will  vary  through  the  bed  due  to 
the  change  in  the  air  temperature  as  it  moves  through  the 
bed.   However,  calculations  based  on  a  change  in  film 
temperature  of  60  F  through  the  bed  would  result  in  a 
variation  in  h  of  only  1%.   Since  this  is  an  upper  limit 
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on  the  variation,  this  effect  can  be  neglected. 

One  might  propose  that  there  could  be  a  variation  in 
h  through  the  bed  caused  by  a  change  in  flow  characteris- 
tics, similar  to  that  which  occurs  for  flow  over  tube 
banks  in  heat  exchangers  (31) .   However,  there  are  no 
indications  of  this  phenomenon  in  any  of  the  literature 
which  was  reviewed  on  heat  transfer  in  deep  beds  of  parti- 
cles.  All  such  studies  used  an  average  h  for  the  entire 
bed  and  in  most  cases  obtained  satisfactory  results. 

Comparison  of  Experimental  Data  With  Numerical  Solution 

The  experimental  results  of  forced-air  cooling  tests 
on  Valencia  oranges  and  Marsh  grapefruit  are  presented  in 
graphical  form  along  with  corresponding  results  from  the 
prediction  model.   Individual  mass  average  and  surface 
temperatures  are  indicated  as  a  function  of  time  and  posi- 
tion within  the  bed  for  each  air  flow  rate.   Appendix  C 
contains  a  complete  set  of  these  comparison  curves,  includ- 
ing all  experimental  data.   Selected  curves  from  this  set 
will  be  included  in  this  section  to  support  discussion. 

The  thermal  and  physical  properties  of  Valencia  oranges 
and  Marsh  grapefruit,  as  listed  in  Table  6,  were  used  to 
determine  input  data  for  the  numerical  solution.   The 
convective  heat  transfer  coefficient  (h)  used  in  the  numer- 
ical solution  was  calculated  from  equation  [84] •   Values  of 
"c"  and  "n"  for  oranges  and  grapefruit  were  taken  from  Table 
6  as  0.374,  0.620,  1.04  and  0.482  respectively. 
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Experimental  data 

The  consistency  of  the  experimental  data,  with  respect 
to  between- fruit  differences  and  between-bulk-load  differ- 
ences, was  considered  to  be  exceptionally  good  when 
compared  to  experimental  temperature  data  collected  by 
other  investigators.   Tests  No.  2  and  No.  6  (Figures  29  and 
30)  were  conducted  on  two  different  bulk  loads  of  oranges 
under  similar  conditions.   Very  little  difference  in  tem- 
perature response  can  be  detected  between  the  tests,  and 
both  sets  of  experimental  data  agree  equally  well  with  the 
theoretical  prediction  curve. 

Each  experimental  data  point  is  an  average  of  two 
values  taken  from  two  fruit  which  were  located  at  the  same 
level  in  the  bed.   An  indication  of  the  between-f ruit 
variation  can  be  obtained  by  noting  the  average  difference 
between  these  values  taken  at  the  approximate  half-cooling 
time.   The  average  difference  was  approximately  two  degrees 
Fahrenheit. 

Analyses  of  the  differences  between  experimental  and 
theoretical  data 

Evaporative  cooling. — The  mass-average  temperatures 
indicated  by  the  prediction  model  fit  the  experimental 
data  reasonably  well.   However,  in  most  cases,  the  differ- 
ence in  temperature  between  the  top  and  bottom  layers  of  the 
bed  was  greater  for  the  experimental  data  than  it  was  for 
the  theoretical  data  (Figure  31) .   The  primary  reason 
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Figure  29.   Comparison  of  experimental  and  theoretical  temp- 
erature response  of  oranges  to  forced-air  cooling 
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Figure  30.   Comparison  of  experimental  and  theoretical  temp- 
erature response  of  oranges  to  forced-air  cooling 
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Figure  31. 


Comparison  of  experimental  and  theoretical  temp- 
erature response  of  oranges  to  forced-air  cooling 
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for  this  is  probably  due  to  a  variation  in  h  through  the 
bed  caused  by  evaporative  cooling  on  the  surface  of  the 
fruit.   Figure  31  indicates  that  the  bottom  two  layers 
cooled  faster  than  predicted  while  the  top  two  layers  cooled 
slower  than  predicted.   This  means  that  the  h  for  the 
bottom  two  layers  was  greater  than  that  for  the  top  two 
layers,  since  the  prediction  curves  were  generated  using 
a  constant  h  through  the  bed.   This  could  be  attributed 
to  evaporative  cooling  since  the  humidity  of  the  air 
increases  as  it  moves  through  the  bed,  causing  less 
evaporative  cooling  in  each  succeeding  layer  of  the  bed. 
An  estimate  of  the  upper  limit  for  the  amount  of  cooling 
produced  by  mass  transfer  can  be  made  by  considering  the 
total  amount  of  moisture  loss  during  cooling.   The  average 
amount  of  moisture  lost  during  a  cooling  test  amounted  to 
0.53%.   If  all  of  the  heat  for  vaporization  came  from  the 
fruit,  the  temperature  of  the  entire  bed  would  be  reduced 
by  approximately  6  F.   However,  a  large  portion  of  the  heat 
comes  from  the  air  and  in  turn  leaves  the  bed  without 
significantly  cooling  the  fruit.   Several  investigators  have 
observed  that,  in  some  cases,  the  temperature  of  the  air 
leaving  the  bed  was  cooler  than  that  entering. 

Radiation  heat  transfer — Heat  was  transferred  from  the 
top  and  bottom  layers  of  the  bed  to  the  cold  inside  sur- 
faces of  the  cooler  by  radiation.   However,  assuming  black 
bodies  and  a  maximum  temperature  difference,  radiation 
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could  account  for  only  about  2  F  of  cooling  during  an 
average  cooling  test.   The  perforated  top  and  bottom  of 
the  box  reduced  the  exposed  surface  area;  thus   the  amount 
of  cooling  by  radiation  would  be  reduced  proportionally. 

Conduction  along  the  thermocouple  wire--Heat  will  be 
conducted  away  from-  the  measuring  junction  in  the  warm 
fruit  because  the  thermocouple  wire  leading  from  the  meas- 
uring junction  is  exposed  to  the  cold  air  surrounding  the 
fruit.   This  will  result  in  a  reading  lower  than  the  true 
value.   Using  basic  heat  transfer  relationships,  the 
dimensions  of  the  wire,  and  average  operating  conductions; 
the  maximum  deviation  was  estimated  to  be  -2  F. 

Insulating  effect  of  rind — The  comparison  curves 
(Appendix  C)  indicate  that  the  experimental  surface 
temperatures  are  consistently  lower  than  the  predicted 
theoretical  values.   Although  some  of  this  difference 
could  be  attributed  to  the  foregoing  reasons,  it  is  believed 
that  the  major  portion  of  it  must  be  attributed  to  some 
other  cause.   Since  the  difference  was  so  much  more  pro- 
nounced on  the  grapefruit  than  on  the  oranges,  the  effect 
of  the  rind  on  the  temperature  gradient  was  further 
investigated.   The  effect  of  the  different  thermal 
properties  of  the  rind  had  previously  been  neglected  based 
upon  experimental  data  in  the  literature. 

In  order  to  investigate  the  effect  of  the  rind  on  the 
temperature  profile  in  citrus  fruit  a  solution  for  heat 
conduction  in  a  composite  sphere  was  needed.   Such  a 
solution  with  convection  boundary  conditions  at  the  surface 
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would  provide  a  comparison  of  surface  temperatures  for 
identical  conditions.   A  solution  with  convection  boundary 
conditions  was  not  found  in  the  literature.   However, 
Carslaw  and  Jaeger  (11)  give  the  following  solution  for  a 
composite  sphere  initially  at  a  temperature,  T1 ,  and  the 
outer  surface  (r=b)  suddenly  lowered  to  zero  temperature 
and  held  constant.   It  is  assumed  that  there  is  no  contact 
resistance  between  layers  (r=a) . 


T.      =  2bT,/r  E  [l/<j>U>]  sin  rX   sin  aX 
inner      1    .    T   n         n       n 
n=l 


"aiXnt 
sin  p  (b-a) Xn  e  [85] 


T      ,         =    2bTn/r   Z    [i/<j>UJ]    sin2aX      sin  p(b-r)X 
outer  r         ,  n  n  ,  n 


-a .  X   t 
e      x   n  [26] 


2  2 

4>  (X    )    =    aaX      sin  pX     (b-a) +p  (b-a)  X      sin   aX 
Yn  n  ^n  n  n 


2  2 

+  [  (l-pa)/paXn]sin   aXn   sin  p(b-a)Xn 


[87] 
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where 


±X  ,  n  =  1,  2,  ...,  are  the  roots  of  equation  [87] 


p  =  /a  .  /a 
^  l'  o 

a  =  k./pk 

a  =  radius  of  inner  core 

b  =  radius  of  outer  core 

T.      =  temperature  of  inner  core 
inner 

T      =  temperature  of  outer  shell 
outer      e 

k.  =  thermal  conductivity  of  inner  core 

k  =  thermal  conductivity  of  outer  shell 
o 

a.  =  thermal  diffusivity  of  inner  core 

a  =  thermal  diffusivity  of  outer  shell 
o 

T,  =  initial  uniform  temperature  of  sphere 
Equations  [85]  through  [87]  were  evaluated  using  a  digi- 
tal computer  and  typical  constants  for  grapefruit. 

Figure  32  gives  a  comparison  of  the  temperature  pro- 
file in  a  homogeneous  sphere  with  average  properties  of 
a  grapefruit  and  the  temperature  profile  in  a  composite 
sphere.   For  the  composite  sphere,  the  inner-core  prop- 
erties were  set  equal  to  those  of  the  juice  vesicles 
and  the  outer-shell  properties  were  set  equal  to  those 
of  the  rind.   These  curves  represent  the  temperature  pro- 
file after  0.5  hours  of  exposure  time.   The  solid  line 
represents  the  temperature  profile  in  the  composite 

sphere  using  a-      =  a   .  _f  but  with  the  thermal  con- 
3   inner    outer' 

ductivity  (k)  differing  by  a  factor  of  two.   The  dashed 


107 


60 


- 


CO 

w 
w 

w 

Q 


D 
E-i 

& 

W 
ft 


Composite  sphere 

Homogeneous  sphere 

Homogeneous  superimposed 

on  composite 
Elapsed  time  -  0.5  hours 


rind 


0.2 


0.4        0.6 
RADIUS  RATIO 


0.8 


1.0 


Figure  32.   Comparison  of  the  temperature  profile  in  a 
homogeneous  sphere  with  average  properties 
of  a  grapefruit  and  the  temperature  profile 
in  a  composite  sphere  with  properties  corre- 
sponding to  the  rind  and  juice  vesicles 
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line  represents  the  temperature  profile  using  the  same 

equation,  but  with  a.      =  a   .    and  k.      =  k   . 

^  inner    outer      inner    outer 

The  temperature  drop  through  the  rind  region  is 
much  greater  for  the  composite  sphere  than  it  is  for  the 
homogeneous  sphere  (without  rind) .   Although  an  actual 
comparison  of  surface  temperatures  cannot  be  made,  an 
estimate  can  be  obtained  by  superimposing  the  profile 
of  one  curve  on  the  other,  starting  at  the  innerface 
between  the  juice  vesicles  and  the  rind.   Thus,  it  is 
seen  that  the  surface  temperature  may  be  as  much  as  10  F 
lower  when  the  effect  of  the  rind  is  considered.   For 
actual  air-cooling  tests  the  effect  of  the  rind  on  the 
temperature  profile  will  not  be  as  great  as  it  was  in  this 
example,  since  the  cooling  rate  will  be  less.   However, 
it  will  no  doubt  account  for  most  of  the  difference  be- 
tween the  experimental  and  predicted  surface  temperatures. 

Suggestions  for  Further  Study 

Optimum  design  and  operating  conditions 

The  results  of  this  study  provide  a  basis  for  making 
an  economic  analysis  of  the  cooling  system.   By  assuming 
certain  cost  figures  and  certain  limitations  based  upon 
the  physiological  response  of  the  product;  the  optimum 
air  flow  rate,  air  temperature,  arrangement  of  product, 
etc. ,  can  be  determined.   The  results  of  this  study  could 
be  in  the  form  of  a  computer  program  so  that  the  designer 
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or  operator  could  plug  in  current  input  data  based  on  his 
particular  application. 

Modification  of  numerical  model 

The  results  of  this  study  are  believed  to  be  suf- 
ficiently accurate  for  engineering  application,  except 
for  the  predicted  surface  temperatures  for  grapefruit. 
Further  work  is  necessary  in  order  to  model  heat  transfer 
in  a  bed  of  composite  spheres. 

In  order  to  apply  the  numerical  model  to  other  pro- 
ducts, it  will  be  necessary  to  account  for  water-vapor 
transfer  from  the  product  and  heat  conduction  in  irregu- 
larly shaped  bodies. 


CONCLUSIONS 

1.  The  numerical  model  can  be  used  to  predict  the 
temperature  distribution  in  an  individual  orange  or 
grapefruit  and  within  the  bulk  load  as  a  function 
of  time. 

2.  For  slow  cooling  rates  (25  fpm  or  less) ,  the  tem- 
perature gradients  within  the  individual  fruit  can 
be  neglected  without  significantly  affecting  the 
predicted  mass-average  temperature.   However,  for 
faster  cooling  rates  the  effect  becomes  significant. 

3.  The  average  convective  heat  transfer  coefficient  (h) 
can  be  predicted  for  a  bed  of  oranges  or  grapefruit 
by  using  the  Nusselt  number  -  Reynolds  number  rela- 
tionship as  indicated  in  equation  [84]  •   The  constants 
in  equation  [84]  can  be  evaluated  by  comparing 
experimental  data  with  data  from  the  numerical  solution, 
thus   providing  a  new  and  direct  method  for  evaluating 

h  in  packed  beds . 

4.  An  average  thermal  diffusivity  (a)  for  oranges  or 
grapefruit  can  be  determined  experimentally  by 
utilizing  a  fast  cooling  method,  such  as  an  agitated 
water  bath.   Since  a  =  pc/k,  all  necessary  properties 
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can  be  determined  with  adequate  accuracy  for  engineer- 
ing calculations  when  a  and  p  are  determined  experi- 
mentally, and  c  is  calculated  as  a  function  of  moisture 
content.   This  allows  k  to  be  calculated,  thus 
eliminating  the  tedious  job  of  experimentally  deter- 
mining k. 

5.  The  presence  of  the  thick  rind  on  grapefruit  signifi- 
cantly affects  the  temperature  gradient  near  the 
surface. 

6.  Radiation  heat  transfer  does  not  significantly  affect 
the  cooling  rate  for  the  normal  range  of  temperature 
used  in  precooling. 
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APPENDIX  A 
COMPUTER  PROGRAM  AND  SAMPLE  OUTPUT  FOR  CONDUCTION 
IN  A  SPHERE  WITH  SPECIFIED  SURROUNDING  TEMPERATURE 

DIMENSION  ATEMP (21)  , PRSNT (21)  ,PREC (21)  ,TESTNO  (3) 
99  READ  (5,1,END=1000)  ALPHA, CONDTY,  HTC , RADIUS ,NPOINT , 
*TPINIT,TIMINC,THETA,INC,  (TESTNO(I)  ,1=1,3) 

1  FORMAT  (F5.5,F5.3,F5.1,F5.3,I5,F5.1,F5.3,F5.2,I5,23X,3A4 

*) 

READ  (5,2)   (ATEMP (L) ,L=1, 11) 

2  FORMAT  (11F5.1) 

PRINT  70,  (TESTNO  (I)  ,1=1,3) 
70  FORMAT  (1H1 ,' UNIVERSITY  OF  FLORIDA  AGRICULTURAL  ENGINE1, 
*'ERING  DEPARTMENT — COOLING  TEST  NUMBER  ',3A4,//) 

PRINT  80,  ALPHA, CONDTY, HTC, RADIUS, NPOINT ,TPINIT,TIMINC , 
*THETA,INC, (ATEMP (L) ,L=1. 11) 
80  FORMAT  (35X, ' THERMAL  DIFFUSIVITY ' , Fll . 5//, 35X , ' THERMAL ' , 
*'  CONDUCTIVITY' ,F8. 3//, 35X, 'CONVECTIVE  COEFFICIENT', 
*F5.1//,35X, 'RADIUS' ,F22. 3//, 35X, 'NUMBER  OF  POINTS ', 19// , 
*F14.3//,35X, 'TOTAL  TIME ' ,F17 . 2// , 35X, ' PRINT  FREQUENCY", 
* 110/////, IX, 'AIR  TEMPERATURES: ' ,11F6.1) 

SEGMNT=NPOINT-l 

RINC=RADIUS/SEGxMNT 

D=l . 0-3 . 0*ALPHA*TIMINC* ( (RADIUS-RINC/2 . 0 ) **2/RIND+RADIU 
*S**2*HTC/CONDTY)/(RADIUS**3- (RADIUS-RINC/2 . 0) ** ) 

IF  (D  .LT.  0.0)  GO  TO  100 

B=3. 0*ALPHA*TIMINC* (RADIUS-RINC/2 . 0 ) **2/RINC/ (RADIUS** 3 
*-  (RADIUS-RINC/2. 0)** 3) 

H=6 . 0*ALPHA*TIMINC/RINC**2 

P=1.0-H 

NDIV=THETA/TIMINC/10 . 0 

TIME=0.0 

NO=0 

N=l 

NUM=1 

TAIR= ATEMP  (1) 

DO  3  J=l, NPOINT 

3  PRSNT (J)=TPINIT 

PRINT  70,  (TESTNO  (I)  ,1=1,3 

PRINT  300 
300  FORMAT  (1H  ,'  TIME   AIR   ****************  PRODUCT  TEMP', 
*'ERATURE  AT  R/RADIUS  **************** ' / ,1X, '   (HR)   TEM ' , 
*'P',31X,  '(F)'/,6X'   (F)    1.0    0.9    0.8    0.7    0.6', 
*'    0.5    0.4    0.3    0.2    0.1    0.0'//) 

PRINT  4,  TIME, TAIR, (PRSNT (J) ,J=1, NPOINT) 

4  FORMAT  (1H  , F6 . 3 ,F5 . 1 , 11F6 . 2 ) 
9  N=N+1 

NCT=0 

DAIR=  (ATEMP  (N)-ATExMP  (N~l)  ) /THETA*  (TIMINC*10  .  0  ) 

5  TIME=TIME+TIMINC 
NCT=NCT+1 
NUM=NUM+1 

DO  15  J=l, NPOINT 
15  PREC (J)=PRSNT (J) 

PRSNT (1)=B*PREC (2 ) +C*TAIR+PREC (1) 
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TAIR=TAIR+DAIR 

M=NP0INT-1 

DO  7  J=2,M 

R=RADIUS-RINC* (J-l) 

E=1.0-ALPHA*TIMINC* (2. 0+ ( (RINC/R) **2/2 . 0 ) )/RINC**2 

IF  (E  .LT.  0.0)  GO  TO  101 

F=ALPHA*TIMINC* (1 . O-RINC/2 . 0/R) **2/RINC**2 

G=ALPKA*TIMINC* (1 . Q+RINC/2 . 0/R) **2/RINC**2 

7  PRSNT (J ) =E*PREC (J ) +F*PRSC ( J+l ) +G*PREC ( J-l ) 
J=NPOINT 

PRSNT (J)=H*PREC (J~1)+P*PREC (J) 
IF  (NUM- INC) 4  0,40,20 
20  PRINT  6,  TIME,  TAIR, (PRSNT (J) , J=l ,NPOINT) 
6  FORMAT  (1H  ,  F6.3,  F5 . 1,11F6 . 2) 

NO=NO+l 
NUM=1 

40  IF  (NO-50)  90,10,10 
10  NO=0 

PRINT  70,  (TESTNO(I)  ,1=1,3) 

PRINT  300 
9  0  IF  (THETA-TIME)  9  9,99,8 

8  IF  (NDIV-NCT)  5,9,5 
100  PRINT  50 

50  FORMAT  (1H1,'D  IS  LESS  THAN  ZERO') 
GO  TO  9  9 

1000  PRINT  1001 

1001  FORMAT  ('1' ) 
STOP 

END 
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APPENDIX  B  -  COMPUTER  PROGRAM  AND  SAMPLE  OUTPUT  FOR  COOLING 

A  DEEP  BED  OF  CITRUS  FRUIT 

01  MTNS  ICN  FP.SNT<21.25),PPEC<21.2M.TESTNC<3),VCL<21>. 

♦  TPSt  CA(21  )  .  TPLSUP(  2  1  )  ,PPfP/>V(  ?  =  )  ,PREPAV<  25) 

qy     P     All      (  S  ,  1   ,FND  =  1  000  )      /LPHA  ,  C'.NOT  Y  ,      I- T  C  ,  P  A  D  I  U  S  .  NPO  I  NT  ,  T  P  I  N 

♦  t  T  .  T  [NINC.ThETA,  INC.  (Tfc"STNC(  I  )  ,1  =  1  ,3) 

1  F 1 1 1. ;  N  A  T     (F5.5.FS.3.F5.2.F5.3.IS.F5.1  .F5.3.FE.2.I5.23X.3A4 

♦  ) 

R'-'AP      (5.2)      ShA,SHP,TCTnT,AIinvTF.NUMAY,BHGHT,ATFVP,AIPVE 

♦L. KPPGC 

2  FGRVAT  (F5.3.F5.3 .F5.1 .F5.0. I5.F5.2 »F5 .1  iF5  .0»  I  1 ) 
PWINT  70.  (  TFSTMCM  I  )  .  1  =  1  .3) 

70  FOPVAT  (lhli'  * . 

♦  •  •  ♦ 3A4 , ////////  ) 
PRINT  «0.   ALFHA  .CONDI Y.HTC  .PAC  IUS  .NPOINT  »TPI  M  T  .  T IM  INC  .T 

+  HrTA.INC.SI-A.SHP.TOThT,AIPteTF.NUN'LAY,F,HGHT,ATEVP,AIRVtrL 
80  POFVAT   (  35X  .  'THfc RMAL  D  I  FF U S  I  V  I  T Y ',  F  1  1  .  5/ /,  35 X ,' THF R MAL '  . 
*•  CONDUCT  IV ITY«  , F6  .3// ,35X  .  'CONVFCT  I  VE  C OEF F  I  C  I  F NT  •  , 
*F6  .  ?//  ,35X  ,  • PAQI US  •  ,F22  .2//.35X  .  • NUVPFf  CF  POINTS'. 
+  IO//.35X,  •  IN  IT  I  AL  FRUIT  TE^P .  •  »F fi .  1  // * 25X  .  * T  I  ME  INCREVE'. 
* 'NT •  ,F 1  4  .3//  ,35X  ,  'TOTAL  T  I  N£  •  .  F  1  7  .  2 //  .  3 5 X ,  •  P P I  NT  FRE •  * 
*«OUENCY«  . I  10//. 35X,«  SPECIFIC  FFAT  A  I P •  ,  F  1  2  .  A// .  3 5 X , 

♦  •SPECIFIC  HEAT  PRODUCT"  .F7.3//.3SX.  • fcE  IGHT  PRODUCT*. 
*F14.  1//.35X,  •  AIR  FLCV\  P  ATE  •  .  F  1  A  .  0//  ,  35  >  .  •  L  A  YERS  IN  BED'. 

♦  I  1 2//. 35X,  »aEO  DEPTH  •  ,F  19.2//  ,25X  .  'ENTERING  AIR  TEMP.'. 
*Fi0.1//.35X,  'AIR  VELCC I T Y •  .  F I  5  .  0/ ) 

IF  (KPPCO-1)  350.350,360 

350  PRINT  351 

351  FORMAT  ( 35X, 'PRODUCT  '. 16X.  'GRAPEFRUIT '  ) 
G  J  T  C  3  5  2 

360  PRINT  361 

361  FORVAT  (  35X  ,  'PPCCUCT ',  16X  , 'ORANGES'  ) 
252  SFJGMNT  =  NPGINT- 1 

PINC-PADIUS/SFGMNT 

D=l  .0-2.0vALFHA*TIMINC*(  (R/>DIUS-RINC/2.0)**2/RINC  +  RA0IUS 
***2*rTC/CCNOTY)/(RADIUSv*3-((<AOIU£-flINC/2.0)**3) 

IF   (C  .LT  .  0.0)  GC  TC  100 

t)=3.0vALPHA*TIMINC*<RA0IUS-RINC/2.0)**2/RINC/(RA0IUS**3- 
#<PA0IUS-RINC/2.0  )  *  *  3  ) 

C  =  3.0*PADIUS**2*I-TC*ALP»-'A*TIMINC/CCNCTY/(RACIUS**3-(RADI 

*US-fi  INC/2 .C )**3) 

b  =  6.0*<ALPt-A*TIMINC/RINC**2 

P=l . 0-H 

NSFG=NPQINT-i 

V,TPL  =  TCTU,T/NUVLAY 

DO  3  1  J=l  ,NSEG 

A  =  J 

Sh'G  =  NSEG 
31     VOL ( J) =<  1  .-( A-l  .  )/ScG  )**3-<  1  .-A/S£G)**2 

IF      (KPPCD-1)     250,250,260 
250     VOLC-VCK1  >*30.0*VCL(2)*30.0 

VOL  1  =  0.0 

CO     200     L=3.NSEG 
20C     VDLI=VCLI+VCL(L)*57.2 

TVCL=VCL  I  +VCL0 

V0L<l)=VCL(l)+(30.0/TV0L) 

121 


122 


20  1 


260 


261 


262 
270 


300 


V""U.  I  «'  )  =VCL  <  2  )  *  (  ,"»0  .0/TVCL  ) 

On  ,:oi  L=3»NSRG 

VUL  (l  >=VCL (L  >*<57.2/TVCL  ) 

C.n   TC  ?7Q 

VOL0  =  VCL  (  1  )  *23  .  7 

VOL  1-0 .  0 

DO  2£1  L=2.NSEG 

VOL I-VCLl+VCL(L)+57.6 

TVCL=VCLI+VCLO 

VOL< 1 )=VCL< 1 )«(23. 7/TVCL) 

CD  262  L-2.NSEG 

VOL  (I )=VCL(L)*(57.6/TV0L) 

T  !  w  fr  =  0  .  0 

N  0  -  0 

N-  1 

TA  IR=ATEMP 

CO  3  J=l  ,NFC  INT 

DO  2  K=  1  ♦NLN'LAY 

PRSNTC  JiK)  =  TPIMT 

PRINT  70,  (TFSTNCH  I  ) ,  1  =  1  .3) 

K=  1 

PFINT  300 

FORMAT  ( 15X. 'DISTANCE   TIME 


AIR        i*********** *****     PRO* 


*»DUCT     TEMPERATURES    AT     R/RACIUS     ***************     NASS*/. 
*15X.'FRCM«,6X,»(t-R)        TEMP •  ,3 1 X t • ( F  )  • • 3 3X • • A VC •/ • 15X « 
*•  INLET*  ,11  X  ,« <F)  1.0  C.9  0.8  C.7  0.6  0.5», 

*•  0.4  0.3  0.2  0.1  0.0        TEMP«/ilSX.  »{FT) • ) 

PRINT     4,     TIME.TA  IS  ,  (FRSNT  <  J  ,  K  }  ,  J=l  ,NPC  INT)  ,TPIMT 

4  FORMAT     (  15X,<;X.Fe.3.F5.  1  .1  1F6.2  ,F6.2> 
CO     37     K=l , NU^LAY 

37     PREPAV(K)=TPIMT 

42  NPRINT=1 

42  NUM=l 

41  TIN'E  =  TIME  +  TIMINC 

K  =  0 

NU«  =  MJ  V  ♦  1 

5  K-K+  1 
A<  =K 

ANU^LY  =  MJN'L  AY 

IP  (K-l  )  2C 1  ,301  ,302 

301  OF  INLT=  (BhC-HT/ANUMLY  >/2  .0 
GO  TG  30  3 

302  DP  INLT-AK*(UHGHT/ANU*'LY)-(eHGHT/ANUMl.Y  J/2.0 

303  00  15  J-l  .NPC I  NT 

IS  PRE C  <  J . K  >  =FR  SNT ( J , K ) 

P;<  S  N  T  (  1  ,K)=H  +  PBl:C(2.K>+C*TAIB-*C  +  PREC(l  ,K) 
C^^CIM-  1 
CO  7  J=2,M 
R=RADIUS-RINC*< J-l) 

E-l   .  G-Al_Pl-A*TIMIiNC*<2.0+<  (RINC/R>**2/2.0))/RlNC**2 
IF  (E  .LT.  0.0)  GO  TC  101 

F=ALPHA*TIKINC*( 1 «0-R INC/2 «0/R ) **2/R INC* *2 
G  =  ALP^A*TIMNC*<1.0+RINC/2.0/h)**2/RINC**2 
7  PRSNT( J,K)=E*P«EC( J,K  )+F*PREC( J+l  ,  K  )  +G *PREC ( J- 1 , K  ) 


123 


J-^ncI^T 

PRSNT*  J.K)=H*PB£C< J- I .K )*P*PREC< J»K ) 

DO  3C  J=l  .NSEG 


30 

32 

33 

20 

6 


DO  3C  J= 1  .NSEG 
TPStCA(J)=(PKSNT<J+l.K)-PRSNT(J,K))*C.5300*PPSNT(J.K> 

CO  3  2  J=l  .NSEG 


TPLStP(J)=1PSECA(J)*VCL(J) 

P  O S  P  A  V ( K  )  =  C.C 
CO  33  J=l  .NSEG 
PPSPAV(K)=PRSPAV(K)-i-TFLSOe<J) 
IF  (NU^-INC)  40.40,20 


IF  (NU^-INC)  40.40.2C 

PRINT  6,DFINLT.TiyE,TAIft,  (FPSNT<J,tO,J=l,NP0INT),PR 
+  K) 


SPAV( 


FORNAT 


AO 
1C 


89 


90 
36 

39 


100 
50 

10  1 
60 

1000 
1001 


(15X,F5.2.4X.F£.3»F5.1.11F6.2.F€.2) 

NPfiIM  =  NPRUT+l 

NG=NC+  1 

IF  (NC-33)  89,10.10 

N0  =  0 

PRINT  70, (TESTNCH  I  )  , 1  =  1  ,3) 

PRINT  300 

DTAIR=WTPL*SHP*(PREPAV(K)-PfiSPAV(K>  >/ < A  I RWTF +T IM INC*SHA ) 

TA IP=TAIR+CTAIR 

IF      (K-NUVLAY)     5.90.90 

IF     (TpETA-TIVE)     99.99,38 

CO     39     K=l  ,  MjVLAY 

PREPAV(K)=FRSPAV(K) 

TAIR^ATEMP 

IF  < NU^LAY-NFPINT  )  43.41,41 

PR  INT  5  0 

FQRVAT  (IhU'O 

GO  TC  9  9 

PRINT  6  C 

FORMAT  (  1H1  *  •£ 

GO  TC  99 

PRINT  100  1 

FORMAT  <•!•) 

STCP 

END 


IS  LESS  THAN  ZERO') 
IS  LESS  THAN  ZEKO  •  ) 
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INPUT  DATA 

THERMAL  DIFFUSIVITY  C.C0437 

THERKAL  C  CtsCLCT  I  V  I  TY  C.24S 
CCNVECTIVE  COEFFICIENT  3.C6C 

RADILS  0.117 

NLMEER  CF  FCINTS  11 

INITIAL  FRLIT  TEMP.  91.0 

TINE  INCREMENT  0.C05 

TOTAL  TIME  2.50 

PRINT  FRECLENCY  20 

SPECIFIC  FEAT  AIR  C.2400 

SPECIFIC  FEAT  PRCDUCT  0.900 

WEIGHT  FRCDLCT  340. C 

AIR  FLCU  PATE  1S3C. 

LAYERS  IN  EEC  11 

BED  DEPTH  2.200 

ENTERING  AIR  TEMP.  32.0 

AIR  VELOCITY  ICC. 

PRGCLCT  ORANGES 
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Figure  33.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  34.   Theoretical  and  experimental  temperature 

response  of  oranges  to  forced-  air  cooling 
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Figure  35.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  36-   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  37.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  38.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  40.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  41.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  42.   Theoretical  and  experimental  termperature 
response  of  oranges  to  forced-air  cooling. 
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Figure  43.   Theoretical  and  experimental  temperature  . 
response  of  oranges  to  forced- air  cooling. 
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Figure  44.   Theoretical  and  experimental  temperature 

response  of  oranges  to  forced-air  cooling. 
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Figure  45.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling. 
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Figure  46.   Theoretical  and  experimental  temperature 

response  of  oranges  to  forced- air  cooling, 
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Figure  48.   Theoretical  and  experimental  temperature 

response  of  oranges  to  forced-air  cooling. 
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Figure  49.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  50.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  51.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  52.   Theoretical  and  experimental  temperature 
response  of  oranges  to  forced-air  cooling 
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Figure  53.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  54.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  55.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  56.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  57.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  58.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  59.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  60. 


Theoretical  and  experimental  temperature 
response  of  grapefruit  to  forced-air  cooling 
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Figure  61.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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Figure  62.   Theoretical  and  experimental  temperature 

response  of  grapefruit  to  forced-air  cooling 
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